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Viscous Incompressible 

Flows 


8.1 INTRODUCTION 

In the analysis of motion of a real fluid, the effect of viscosity should be given 
consideration. Influence of viscosity is more pronounced near the boundary of a 
solid body immersed in a fluid in motion. The relationship between stress and rate 
of strain for the motion of real fluid flow was first put forward by Sir Isaac 
Newton and for this reason the viscosity law bears his name. Later on, G.G. 
Stokes, an English mathematician and C.L.M.H. Navier, a French engineer, 
derived the exact equations that govern the motion of real fluids. These equations 
are in general valid for compressible or incompressible laminar flows and known 
as Navier-Stokes equations. When a motion becomes turbulent, these equations 
are generally not able to provide with a complete solution. Usually, in order to 
obtain accurate results for such situations, the Navier-Stokes equations are 
modified and solved based on several semi-empirical theories. In the recent past, 
some researchers have proposed that, on a fine enough scale, all turbulent flows 
obey the Navier-Stokes equation and computationally, if a fine enough grid is 
used with appropriate discretization methods, may be both the fine scale and large 
scale aspects of turbulence can be captured. However, in this chapter we shall 
discuss the equation of motion for laminar flows and various other aspects of 
laminar incompressible flows. 
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8.2 GENERAL VISCOSITY LAW 

The well-known Newton’s viscosity law is 

X= ^~L M 

where n is the coordinate direction normal to the solid-fluid interface, ,U is the 
coefficient of viscosity and V is velocity. This law is valid for parallel flows. 
There are more generalized relations which can relate stress field and velocity 
field for any kind of flow. Such relations are called constitutive equations. We 
shall consider here the Stokes’ viscosity law. 

According to Stokes’ law of viscosity, shear stress is proportional to rate of shear 
strain so that 
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The first subscript of T denotes the direction of the normal to the plane on 
which the stress acts, while the second subscript denotes direction of the force 
which causes the stress. 

The expressions of Stokes’ law of viscosity for normal stresses are 
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where /l' is related to the second coefficient of viscosity jl ] by the relationship 
2 

fl' = - — (jl — id ] ). We have already seen that thermodynamic pressure 
p=- -——-—-—. Now, if we add the three Eqs (8.3a), (8.3b) and (8.3c), 


we obtain 
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For incompressible fluid, = V- V = 0 

dx dy dz 


So ,p = - 




yy 


° zz ) . 


is satisfied in the same manner. For compressible 


fluids, Stokes’ hypothesis is // = - — ,Ll. Invoking this to Eq. (8.4), will finally 


conclude that p = 




yy 


<* zz ) 


. Generally, fluids obeying the ideal gas 


equation follow this hypothesis and they are called Stokesian fluids. It may also 
be mentioned that the second coefficient of viscosity, has been verified to be 
negligibly small. 

Now, we can write 
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(8.5a) 


(8.5b) 


(8.5c) 


In deriving the above stress-strain rate relationship, it was assumed that a fluid 
has the following properties 

1. Fluid is homogeneous and isotropic, i.e. the relation between 
components of stress and those of rate of strain is the same in all 
directions. 

2. Stress is a linear function of strain rate. 

3. The stress-strain relationship will hold good irrespective of the 
orientation of the reference coordinate system. 

4. The stress components must reduce to the hydrostatic pressure 
(typically thermodynamic pressure = hydrostatic pressure) p when all 
the gradients of velocities are zero. 


8.3 NAVIER-STOKES EQUATIONS 


Generalized equations of motion of a real flow are named after the inventors of 
them and they are known as Navier-Stokes equations. However, they are derived 
from the Newton’s second law which states that the product of mass and 
acceleration is equal to sum of the external forces acting on a body. External 
forces are of two kinds—one acts throughout the mass of the body and another 
acts on the boundary. 

The first one is known as body force (gravitational force, electromagnetic 
force) and the second one is surface force (pressure and frictional force). 
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Let the body force per unit mass be 


ft = i f x + j f y + kf z 


( 8 . 6 ) 


and surface force per unit volume be 


F = i F x + j F y +k F, 


(8.7) 


Consider a differential fluid element in the flow field (Fig. 8.1). We wish to 
evaluate the surface forces acting on the boundary of this rectangular parallel¬ 
epiped. 



x 

Fig. 8.1 Definition of the components of stress and their locations in a 
differential fluid element 

To accomplish this, we shall consider surface force on the surface AEFfD, per 
unit area, 



Surface force on the surface BFGC per unit area is 


Net force on the body due to imbalance of surface forces on the above two 
surfaces is 


dF sx 


dr d y dz 


dx 


( 8 . 8 ) 
























The McGraw-Hill Companies 


| 318 


Introduction to Fluid Mechanics and Fluid Machines 


Total force on the body due to net surface forces on all six surfaces is 


dF dF. 

VS -l r. r $ 


sy 


dF- 


dx dy dz 
and the resultant surface force dF per unit volume, is 


dr dj’ dz 


(8.9) 
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( 8 . 10 ) 


dx dy dz 

The quantities F sx , F sv and F s: are vectors which can be resolved into normal 
stresses denoted by a and shearing stresses denoted by x as 
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( 8 . 11 ) 


The stress system is having nine scalar quantities. These nine quantities form a 
stress tensor. The set of nine components of stress tensor can be described as 
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( 8 . 12 ) 


The above stress tensor is symmetric, which means that two shearing stresses 
with subscripts which differ only in their sequence are equal. Considering the 
equation of motion for instantaneous rotation of the fluid element (Fig. 8.1) 
about y axis, we can write 

ffl V df y = ( T *_- d >’ dz ) d - Y - ( T ZX dv d y)fc 
= (Cz - F) dV 

where dV is the volume of the element, and 0) v and d/ are the angular accel¬ 
eration and moment of inertia of the element abouty-axis respectively. Since d I 
is proportional to fifth power of the linear dimensions and dV is proportional to 
the third power of the linear dimensions, the left hand side of the above equation 
vanishes faster than the right hand side on contracting the element to a point. 
Flence, the result is 

Fz ^zx 

From the similar considerations about other two remaining axes, we can 
write 

Fy~ Fx 
Fz Fy 

which has already been observed in Eqs (8.2a), (8.2b) and (8.2c) earlier. 
Invoking these conditions into Eq. (8.12), the stress tensor becomes 
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(8.13) 


Combining Eqs (8.10), (8.11) and (8.13), the resultant surface force per unit 
volume becomes 
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As per the velocity field. 
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(8.15) 


By Newton’s law of motion applied to the differential element, we can write 
p(dx dy dz) = (clF) (dx dy dz) + pf h (dr dy dz) 


D t 


or 




Substituting Eqs (8.15), (8.14) and (8.6) into the above expression, we obtain 
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(8.16a) 

(8.16b) 

(8.16c) 


In order to express anc | i n terms of field derivatives, Eqs (8.2) 

Dt Dt Dt 

and (8.5) are introduced into Eq. (8.16) and we obtain 
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These differential equations are known as Navier -Stokes equations. At this 
juncture, it is necessary to discuss the equation of continuity as well, which is 
having a general form 


dp + d(pu) | d(pv) | d(pw) =Q 
dt dx dy dz 


(8.18) 


The general form of continuity equation is simplified in case of incompressible 
flow where p =constant. Equation of continuity for incompressible flow becomes 
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Invoking Eq. (8.19) into Eqs (8.17a), (8.17b) and (8.17c), we get 
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In short, vector notation may be used to write Navier-Stokes and continuity 
equations for incompressible flow as 

DV 
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= P/ft -V P + pV-V 


( 8 . 21 ) 


and 
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( 8 . 22 ) 






















































The McGraw-Hill Companies 



Viscous Incompressible Flows 


OH 


We observe that we have four unknown quantities, u, v, w and p, and four 
equations,—equations of motion in three directions and the continuity equation. 
In principle, these equations are solvable but to date generalized solution is not 
available due to the complex nature of the set of these equations. The highest 
order terms, which come from the viscous forces, are linear and of second order. 
The first order convective terms are non-linear and hence, the set is termed as 
quasi-linear. 

Navicr Stokes equations in cylindrical coordinate (Fig. 8.2) are useful in 
solving many problems. If v n v 0 and v : denote the velocity components along the 
radial, cross-radial and axial directions respectively, then for the case of 
incompressible flow, Eqs (8.21) and (8.22) lead to the following system of 
equations: 



Fig. 8.2 Cyl + indrical polar coordinate and the velocity components 
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(8.23a) 


(8.23b) 


(8.23c) 

(8.24) 


Let us quickly look at a little more general way of deriving the Navier-Stokes 
equations from the basic laws of physics. Consider a general flow field as 
represented in Fig. 8.3. Let us imagine a closed control volume V 0 within the 






























The McGraw-Hill Companies 


322 | | Introduction to Fluid Mechanics and Fluid Machines 


flow field. A control surface, A 0 , is defined as the surface which bounds the 
volume V 0 . The control volume is fixed in space and the fluid is moving through 
it. The control volume occupies reasonably large finite region of the flow field. 

According to Reynolds transport theorem, we know that the laws of physics 
which are basically stated for a system, can be re-stated for a control volume 
through some integral relationship. However, for momentum conservation, the 
Reynolds transport theorem states, “The rate of change of momentum for a system 
equals the sum of the rate of change of momentum inside the control volume and 
the rate of efflux of momentum across the control surface.” 

Again, the rate of change of momentum for a system (in our case the control 
volume is the system) is equal to the net external force acting on it. Now, we shall 
transform these statements into equation by accounting for each term. 



Fig. 8.3 Finite control volume fixed in space with the fluid moving through it 


Rate of change of momentum inside the control volume 

: fffpPdV 


= d_ 
dt. 


Ill dV (since t is independent 


Vo of space variable) 

Rate of efflux of momentum through control surface 

= jjpF(f-dA) = jjpFF-ndA 
A A 


(8.25) 
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IJci-di [dis symmetric stress tensor] 


= jj/(V- CT )dV 

Vo 

and, the body force acting on the control volume 

~l\\pf dv 

V 0 

/ in Eq. (8.28) is the body force per unit mass. 
Finally, we get, 
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We know that —- + V-pV = 0 is the general form of mass conservation 
dt 

equation, valid for both compressible and incompressible flows. Invoking this 
relationship in Eq. (8.29), we obtain 
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Equation (8.30) is often referred to as Cauchy’s equation of motion. In this 
equation, < 7 , the stress tensor, is given by 


(7 = —pi* + p! (V-F) + 2p (Def V ) 


(8.31) 
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From Eq. (8.31), we get 

V<7= -Vp + in' + n) V (V-f) +/tV 2 F (8.32) 


2 

and also from Stokes’s hypothesis, /u' H- fJL = 0 (8.33) 

Invoking Eq. (8.32) in Eq. (8.30) and introducing Eq. (8.33) will yield 


p^- =-Vp + /uV 2 V + ^iuV(V-V) + pf 


(8.34) 


This is the most general form of Navier-Stolces equation. The specific forms 
for different coordinate systems can easily be obtained from Eq. (8.34). 

In a cartesian coordinate system, 
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This is known as deformation tensor. It can be readily seen that (Def V) is a 
symmetric tensor. 


8.4 EXACT SOLUTIONS OF NAVIER-STOKES 
EQUATIONS 


The basic difficulty in solving Navier-Stokes equations arises due to the 
presence of nonlinear (quadratic) inertia terms on the left hand side. However, 
there are some nontrivial solutions of the Navier-Stokes equations in which the 
nonlinear inertia terms are identically zero. One such class of flows is termed as 
parallel flows in which only one velocity term is nontrivial and all the fluid 
particles move in one direction only. 

Let us choose x to be the direction along which all fluid particles travel, i.e. 
u 0, v = w = 0. Invoking this in continuity equation, we get 


du_ + dv^_ + dvf_ =0 
dx dfy pz 



dx 


which means u = it (y, z, t ) 


Now, Navier-Stokes equations for incompressible flow become 
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So, we obtain 


and 


dp dp 

= —— = 0, which means p = p(x) alone, 
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(8.35) 


dy 2 dz 2 

For a steady two dimensional flow through parallel plates, Eq. (8.35) is further 
simplified and analytical solution can be obtained. An insightful description of 
many such analytical solutions in different geometrical configurations has been 
well documented in White [1] and Faber [2]. Flowever, we shall discuss some of 
the important exact solutions in the following sections. 


8.4.1 Parallel Flow in a Straight Channel 

Consider steady flow between two infinitely broad parallel plates as shown in 
Fig. 8.4. Flow is independent of any variation inz direction, hence, zdependence 
is gotten rid of and Eq. (8.35) becomes 



Fig. 8.4 Parallel flow in a straight channel 
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The boundary conditions are at y = b,u = 0; and y 
From Eq. (8.36), we can write 
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(8.36) 
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Applying the boundary conditions, the constants are evaluated as: 


C, = 0 


and 


c = _1 dp b~ 

jj. dr 2 


So, the solution is 


1 d/? 2 2 \ 

u = --— — (b -y ) 
2/u dr 


(8.37) 


which implies that the velocity profile is parabolic. We can extend our analysis 
little further in order to establish the relationship between the maximum velocity 
and average velocity in the channel. 

At y = 0, u = t/ max ; this yields 
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(8.38c) 


The shearing stress at the wall for the parallel flow in a channel can be 
determined from the velocity gradient as 


f n.A 


y*\b 


du 
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-b—--2n 
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Since the upper plate is a “minus y surface”, a negative stress acts in the positive 
x direction, i.e. to the right. 

The local friction coefficient, Ch is defined by 
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(8.38d) 
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where Re = U, lv (2h)iv is the Reynolds number of flow based on average velocity 
and the channel height (2b). Experiments show that Eq. (8.38d) is valid in the 
laminar regime of the channel flow. The maximum Reynolds number value 
corresponding to fully developed laminar flow, for which a stable motion will 
persist, is2300.Inareasonably careful experiment, laminar flow can be observed 
up to even Re = 10,000. But the value below which the flow will always remain 
laminar, i.e. the critical value of Re is 2300. 


8.4.2 Couette Flow 


Another simple solution for Eq. (8.35) is obtained for Couette flow between two 
parallel plates (Fig. 8.5). Elere, one plate is at rest and the other is moving with a 
velocity U. Let us assume the plates are infinitely large in z direction, so the 


z dependence is not there and the governing equation is —— = 

dr 

to the boundary conditions at y = 0, u = 0 and y = h, u = U. 

U 


d 2 u 

dy 2 


subjected 



Fig. 8.5 Couette flow between two parallel flat plates 


We get, 


1 dp 2 , n 

——y +C x y + C 2 

2 p dr 


Invoking the condition (at y = 0, u = 0), C 2 becomes equal to zero. 

1 dp 2 . n 

u = - y + C,y 

2p dx 1 

Invoking the other condition (at y=h,u= U), 
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1 h 2 p dr 


2p dr 

Equation (8.39) can also be expressed in the form 


(8.39) 
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or 


where 
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— = - + P- 1-- 
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u 

p = - 


h h 
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2jj.U 


dp 

V dxy 


(8.40a) 


Equation (8.40a) describes the velocity distribution in non-dimensional form 
across the channel with P as a parameter known as the non-dimensional pressure 
gradient. When P = 0, the velocity distribution across the channel is reduced to 


u_ „ y_ 

U h 

This particular case is known as simple Couette flow. When P > 0, i.e. for a 
negative or favourable pressure gradient (- dp/dx) in the direction of motion, the 
velocity is positive over the whole gap between the channel walls. For negative 
value ofP(P<0), there is a positive or adverse pressure gradient in the direction 
of motion and the velocity over a portion of channel width can become negative 
and back flow may occur near the wall which is at rest. Figure 8.6a shows the 
effect of dragging action of the upper plate exerted on the fluid particles in the 
channel for different values of pressure gradient. 
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Fig. 8.6a Velocity profile for the Couette flow for various values of 
pressure gradient 


Maximum and Minimum Velocities The quantitative description of non- 
dimensional velocity distribution across the channel, depicted by Eq. (8.40a), is 
shown in Fig. 8.6b. The location of maximum or minimum velocity in the channel 
is found out by setting the derivative du/dy equal to zero. From Eq. (8.40a), we 
can write 
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h h 




For maximum or minimum velocity, 
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which gives 


d u _ o 
dy 

Z = 1 

h 2 


1 

2 P 


(8.40b) 


It is interesting to note that maximum velocity for/ 1 = 1 occurs at y/h = 1 and 

equals to U. For P > 1, the maximum velocity occurs at a location — < 1. This 

h 

means that with P> 1, the fluid particles attain a velocity higher than that of the 

moving plate at a location somewhere belowthemovingplate. Ontheother hand, 

when/’ = —1, the minimum velocity occurs, according to Eq. (8.40b), at — = 0. 

y ^ 

For P < —1, the minimum velocity occurs at a location — >0. This means that 

h 

there occurs a back flow near the fixed plate. The values of maximum and 
minimum velocities can be determined by substituting the value of y from 
Eq. (8.40b) into Eq. (8.40a) as 


Hrrt in 


f/(l + P) 2 
4 P 

U(l + P) 2 
4 P 


for P > 1 


for P < 1 


(8.40c) 



-0.2 .0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 


u 

U 

Fig. 8.6b Velocity distribution of the Couette flow 

8.4.3 Hagen Poiseuille Flow 

Consider fully developed laminar flow through a straight tube of circular cross- 
section as in Fig. 8.7. Rotational symmetry is considered to make the flow two- 
dimensional axisymmetric. Let us takez-axis as the axis of the tube along which 
all the fluid particles travel, i.e. 

v : * 0, v r = 0, v 9 = 0 
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Fig. 8.7 Hagen-Poiseuille flow through a pipe 

Now, from continuity equation, we obtain 
vjf + dv, 

r dz 
dv 


dv° 

dr 


= 0 


For rotational symmetry, — • = 0 

r dO 


dz 


Invoking 


= 0 which means v, = v. (r, t) 


dv 


v r - 0, v () - 0, —— = 0, and — (any quantity) = 0 


dz dO 

Navier-Stokes equations, we finally obtain 


in the 


dv z 

dt 


J_ dp 
p dz 


d V z 
dr 2 


1 dv z 

r dr 


For steady flow, the governing equation becomes 


d~v, 1 dv. 


p dz 


dr r dr 
The boundary conditions are 

at r = 0, v. is finite and 
at r = R,v z = 0. 

Equation (8.42) can be written as 


(8.41) 


(8.42) 


or 


or 


or 


or 


d 2 v z 

dr 2 

_d 

dr 


dv z _ 1 dp 


dr 

dv z 

dr 

du_ 
r — 
dr 

dv : 

dr 


v. = 


p dz 

1 dp 

-— i 

p dz 
1 dp 

2 p dz 

1 dp 
2 p dz 

1 dp 
4 p dz 


— r + A 


r 

+ A In r + B 
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at r = 0, v z is finite which means A should be equal to zero and at r = R, v z = 0 
yields 

B =- — ^-R 2 


R z 


v- = 


4 p dz 

dp 


4/uy dz 


v R-j 


(8.43) 


This shows that the axial velocity profile in a fully developed laminar pipe 
flow is having parabolic variation along r. 

At r = 0, as such, v. = v. 


max 

-max 4^^ fc 


\ 

The average velocity in the channel, 

R 


(8.44a) 


v. = 


Q_ 

nR 2 


^2nrv z {r) dr 


nR~ 


or 


or 


v. = 


2jr— ( 

dp) 

R 2 R 4 

4 ju 

dz) 

1 

N 

(N 

_1 


nR 1 


R 1 


dp 


'• ‘i ai-trr- 

v, =2v. 


Now, the discharge through a pipe is given by 


Q = nR- v z 


or 


or 


Q = nR~ 


R 1 


dp 
8p V dz 


[From Eq. 8.44b] 


c _ nP 4 (d P \ 

J 128 pldz, 


(8.44b) 

(8.44c) 

(8.45) 


(8.46) 


Equation (8.46) is commonly used in the measurement of viscosity with the help 
of capillary tube viscometers. Such a viscometer consists of a constant head tank 
to supply liquid to a capillary tube (Fig. 8.8). 


C - Capillary 

M - Manometer 

N - Measuring cylinder 

P - Pump 

R - Reservoir 

T - Tank (constant head) 

V- Valve 



Fig. 8.8 Schematic diagram of the experimental facility for determination of viscosity 
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Pressure drop readings across a specified length in the developed regionof the 
flow are taken with the help of a manometer. The developed flow region is 
ensured by providing the necessary and sufficient entry length. From Eq. (8.46), 
the expression for viscosity can be written as 


nD 4 dp 77D 4 (pi~p 2 ) 


M 128 0 d/ 128 0 l 

The volumetric flow rates (0) are measured by collecting the liquid in a 
measuring cylinder. The diameter ( D ) of the capillary tube is known beforehand. 
Now the viscosity of a flowing fluid can easily be evaluated. 

Shear stress profile across the cross-section can also be determined from this 
information. Shear stress at any point of the pipe flow given by 

dv 


T I, = P 


dv- 


dr 

2 -h 2 


R~ (dp 


From Eq. (8.43), —— = - — 


or 


which means 


dr 

dtn 

dr 


4/i ^ dz 
1 f dp 




R 


U. = - 


2p V dz 
1 f dp 
dz 


(8.47a) 

(8.47b) 


This also indicates that r varies linearly with the radial distance from the axis. 
At the wall, r assumes the maximum value. 


At r = R, 


1 


dp 

ydzj 


R 


Again, over a pipe length of l, the total shear force is 


F , = W 2 jtR-I 


or 


F s = - 
s 2 


Pi~ Pi 

l 


2n R 2 -l 


y 


or F s = -JtR x [Pressure drop between the specified length] 

as it should be. Negative sign indicates that the force is acting in opposite to the 
flow direction. 

Flowever, from Eq. (8.44b), we can write 


1 (dp 




or 


8^(^_-)a 

R 2 


Over a finite length /, the head loss h f = 


pressure drop 
PS 


(8.47c) 

(8.48) 

(8.49) 
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Combining Eqs (8.48) and (8.49), we get 


or 


h f = 


h f = 


32 p{v : ) 


1 


D- 


v /■ — 
Pg 


32 p{v z fl ! 


D pg ( 


(8.50) 


On the other hand, the head loss in a pipe flow is given by Darcy-Weisbach 
formula as 


h f = 


2gD 


(8.51) 


where “/’ is Darcy friction factor. Equations (8.50) and (8.51) will yield 

32 p(v z )lj 1 _ fl(v z )l % , 


D~Pg (P_-)a 


2 gD 


which finally gives /= , where Re = av ^ is the Reynolds number. 

Re n 

So, for a fully developed laminar flow, the Darcy (or Moody) friction factor is 
given by 


64 _ 
' Re 


(8.52a) 


Alternatively, the skin friction coefficient for Hagen-Poiseuille flow can be 
expressed by 


r _ Patr = R\ 
f yP(R-)av 

With the help of Eqs (8.47b) and (8.47c), it can be written 

C f = — (8.52b) 

J Re 

The skin friction coefficient C is called as Fanning’s friction factor. From 
comparison of Eqs (8.52a) and (8.52b), it appears 

f = 4C f 

For fully developed turbulent flow, the analysis is much more complicated, 
and we generally depend on experimental results. Friction factor for a wide range 
of Reynolds number (10 4 to 10 8 ) can be obtained from a look-up chart which we 
shall discuss later. Friction factor, for high Reynolds number flows, is also a 
function of tube surface condition. However, in circular tube, flow is laminar for 
Re < 2300 and turbulent regime starts with Re > 4000. In between, transition 
from laminar to turbulent is induced. 
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As it has been pointed out, the surface condition of the tube is another 
responsible parameter in determination of friction factor. Friction factor in the 


turbulent regime is determined for different degree of surface-roughness 


V A 

of the pipe, where £ is the dimensional roughness and D h is usually the hydraulic 
diameter of the pipe. Friction factors for different Reynolds number and surface- 
roughness have been determined experimentally by various investigators and the 
comprehensive results are expressed through a graphical presentation which is 
known as Moody Chart after L.F. Moody who compiled it. This will be presented 
in detail in Chapter 11. 

The hydraulic diameter which is used as the characteristic length in 
determination of friction factor, instead of ordinary geometrical diameter, is 
defined as 


A =^jr~ (8-53) 

where A w is the flow area and P w is the wetted perimeter. 

Now, let us quickly look at the so called kinetic energy correction factor (a) 
and momentum correction factor (/}) associated with fully developed laminar 
flow. 


Kinetic energy correction factor, a The kinetic energy associated with the fluid 


flowing with its profile through elemental area cL4 = 


-p^-A A 


ZC and the total 


kinetic energy passing through per unit time 


P 

2 



dA. This can be related to 


the kinetic energy due to average velocity (u.) av , through a correction factor, a, as 


ArX A 


p r 3 

a = — 1 ” 


or 


a= — 
A 


2 \v\ dA 

ill 




Flere, for Flagen-Poiseuille flow, 


a= — 
A 


-J 

A J 


Oz) n 


cL4 


1- 


R 


2 A 


Oz)r 


In rdr = 2 


(8.54a) 


Momentum correction factor, p The momentum associated with fluid flowing 
with its profile through elemental area cL4 = p y, 2 cL4; and the total momentum 

passing through any particular section per unit time = p \vl cL4. This can be 
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related to the momentum rate due to average flow velocity (u,) av . through a 
correction factor /j, as 

l2 

dA 


[p(v : ) 2 av A]p = p jv; dA or P ~ 
Here, for Hagen-Poiseuille flow, 


(w)a 


1 R 

?=-J 




1- 


2 A 


Oz)r 


2n rdr = — 
3 


(8.54b) 


8.4.4 Flow between Two Concentric Rotating Cylinders 

Another example which leads to an exact solution of Navier-Stokes equation is 
the flow between two concentric rotating cylinders. Consider flow in the annulus 
of two cylinders (Fig. 8.9), where r x and r 2 are the radii of inner and outer 
cylinders, respectively, and the cylinders move with different rotational speeds 
ft)j and co 2 , respectively. 



'•i 


r 2 



Fig. 8.9 Flow beween two concentric rotating cylinders 

From the physics of the problem we know, v : = 0, v r = 0. From the continuity 
Eq. (8.24) and these two conditions, we obtain 
dv e 


de 


= o 


which means v g is not a function of 9. We assume zdimensionto be large enough 

so that end effects can be neglected and (any property) = 0. Now, we can say 

dz 

v 0 = v 0 (r). With these simplifications and assuming that “0 symmetry” holds 
good, (8.23) reduces to 


vjt = dp_ 
r dr 


and 


d 2 v, 


e + l. dv e 




(8.55) 

(8.56) 


dr 2 r dr 

Equation (8.55) signifies that the centrifugal force is supplied by the radial 
pressure, exerted by the wall of the enclosure on the fluid. In other words, it 
describes the radial pressure distribution. 
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From Eq. (8.56), we get 


_d 

dr 


1 d, ; 

-—(rv e ) 
r dr 


or 


— (rv q )=Ar or v R = 
dr V 9 9 2 


= 0 

Ar , B 


(8.57) 


For the azimuthal component of velocity, Vg, the boundary conditions are: at 
r = r,, Vq = /•) 0 ) ] at r = r 2 , v e = r 2 0 ) 2 . Application of these boundary conditions on 
Eq. (8.57) will produce 


and 


zl= 2 


B = 


co i 


1 07 ->i 2 ) 


(®i -co 2 ) 


2 2 
>\ r 2 


(CO 1 -co 2 ) 


(ri~r{) 

Finally, the velocity distribution is given by 


1 


v a = 


6 (A-'f) 


2 „ 2\ , r \ r 2 (CO\ ~C0 2 ) 

(ft), r 2 — (Wj i\ ) r H- 


Now, T ,. g = B Y rO i s general stress-strain relation. 


or 


In our case. 


or 


\e = B 
*re = B 


f dv e 1 dv r v e 


dr r 39 

_ 

dr r 


\re = B '■ 


V w ' ' J 

Cl ( Vg 


dry r 

Equations (8.58) and (8.59) yields 

2/u 


\-o = 


(co 2 - co x )r 1 2 r 2 \ 


(8.58) 


(8.59) 


(8.60) 


Now, 


2 jiri (co 2 


1 rO\r=r 


and 


2Bn (co 2 -®i) 


l r0\ 


For the case, when the inner cylinder is at rest and the outer cylinder rotates, 
the torque transmitted by the outer cylinder to the fluid is 


_ 2 bK co 

L 2 


2 - Y^'2nr 2 lr 2 


or 


T 2 = An bI 


2 2 

n r 2 

r-f - rf 


2®2 


(8.61) 
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where / is the length of the cylinder. The moment 7\, with which the fluid acts on 
the inner cylinder has the same magnitude. If the angular velocity of the external 
cylinder and the moment acting on the inner cylinder are measured, the coefficient 
of viscosity can be evaluated by making use of the Eq. (8.61). 

8.5 LOW REYNOLDS NUMBER FLOW 

We have seen in Chapter 6 that Reynolds number is the ratio of inertia force to 
viscous force. For flow at low Reynolds number, the inertia terms in the 
Navier-Stokes equations become small as compared to viscous terms. As such, 
when the inertia terms are omitted from the equations of motion, the analyses are 
valid for only Re « 1. Consequently, this approximation, linearizes the Navier- 
Stokes equations and for some problems, makes it amenable to analytical 
solutions. We shall discuss such flows in this section. Motions at very low 
Reynolds number are sometimes referred to as creeping motion. 

8.5.1 Theory of Hydrodynamic Lubrication 

Thin film of oil, confined between the interspace of moving parts, may acquire 
high pressures up to 100 MPa which is capable of supporting load and reducing 
friction. The salient features of this type of motion can be understood from a study 
of slipper bearing (Fig. 8.10). The slipper moves with a constant velocity U past 
the bearing plate. This slipper face and the bearing plate are not parallel but 
slightly inclined at an angle of a. A typical bearing has a gap width of 0.025 mm 
or less, and the convergence between the walls may be of the order of 1/5000. It is 
assumed that the sliding surfaces are very large in transverse direction so that the 
problem can be considered two-dimensional. 



For the analysis, we may assume that the slipper is at rest and the plate is 
forced to move with a constant velocity U. The height h(x) of the wedge between 
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the block and the guide is assumed to be very small as compared with the length / 
of the block. This motion is different from that we have considered while 
discussing Couette flow. The essential difference lies in the fact that here the two 
walls are inclined at an angle to each other. Due to the gradual reduction of 

narrowing passage, the convective acceleration u - is distinctly not zero. 

dx 

However, a relative estimation of inertia term with respect to viscous term 
suggests that, for all practical purposes, inertia terms can be neglected. The 
estimate is done in the following way 

Inertia force _ pu(du/dx) _ pU 2 /I _pUl( h Y 
Viscous force p(d 2 u/dy 2 ) pU/h 2 p W, 

The inertia force can be neglected with respect to viscous force if the modified 
Reynolds number, 


. = mfh f 

v It, 


«i 


The equation for motion iny direction can be omitted since the v component of 
velocity is very small with respect to u. Besides, in the x-momentum equation, 
d 2 u/dx 2 can be neglected as compared with d 2 u/dy 2 because the former is 
smaller than the latter by a factor of the order of (h/I) 2 . With these simplifications 
the equations of motion reduce to 


d 2 u d p 

dy 2 dx 

The equation of continuity can be written as 

h(x) 

Q= \ u dy 
o 

The boundary conditions are: 
at y = 0, u = U at x = 0, p = p 0 

at y = h, u = 0 and at x = I, p = p 0 


(8.62) 


(8.63) 


(8.64) 


Integrating Eq. (8.62) with respect to y, we obtain 

Q 


1 dp 2 

u ~-——y + Qy ■ 

2 p dx 


Application of the kinematic boundary conditions (at y ~ 0, u — U and at y = 
u = 0), yields 


u = U\ \ — 


dp \ , 


h dp [ l y 
2 p dx 


y 


(8.65) 


It is to be noticed that — is constant as far as integration along y is concerned, 
\dx) 

but p and —^— vary along x-axis. At the point of maximum pressure, —— = 0, 
dx dx 

hence 
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U=U I 1 - - 


( 8 . 66 ) 


Equation (8.66) depicts that the velocity profile alongy is linear at the location 
of maximum pressure. The gap at this location may be denoted as h . 

Now, substituting Eq. (8.65) into Eq. (8.63) and integrating, we get 


Q= 

2 


p'h 3 

12 fi 


or 




(8.67) 


where p' = dp/dx. 

Integrating Eq. (8.67) with respect to x, we obtain 


Jfo* - 


dx 


(/?! — ax) 

6 pU 6 pQ 


\2pQ\ 


J * ,+c, 

J (7?! - ax) 3 

(8.68a) 

- + Q 

(8.68b) 


or p = 

a{h x -ax) a(h x -ax) 
where a= (h x — h 2 )/l and C 3 is a constant. 

Since the pressure must be the same (p =p 0 ), at the ends of the bearing, namely, 
p = p 0 at x = 0 and p = p 0 at x = /, the unknowns in the above equations can be 
determined by applying the pressure boundary conditions. We obtain 

Uh ih 2 , 6 fiU 


Q = 


h\ + h 2 


and 


C 3 Po ~ 


a{h x +h 2 ) 


With these values inserted, the equation for pressure distribution (8.68)becomes 

6 pUx (h - h 2 ) 


P = P 0 + 


or 


P~P 0 


h z {h x +h 2 ) 
6pUx (h - h 2 ) 


h 2 (h x 


h 2 ) 


(8.69) 


It may be seen from Eq. (8.69) that, if the gap is uniform, i.e. h = h x = h 2 , the 
gauge pressure will be zero. Furthermore, it can be said that very high pressure 
can be developed by keeping the film thickness very small. Figure 8.10shows the 
distribution of pressure throughout the bearing. 

The total load bearing capacity per unit width is 




6 p, U r x{h — h 2 ) 


2 0 


Zf 


dx 


After substituting h = h x - ax with a = (h , - h 2 )/l in the above equation and 
performing the integration, 


P = 


6pUI 2 
(h x — h 2 )~ 


In — ■ 


h\ ~ f h 
A 1 +fl 2. 


(8.70) 


The shear stress at the bearing plate is 





















The McGraw-Hill Companies 


340 


Introduction to Fluid Mechanics and Fluid Machines 


Tn = 


du 

% 


P - + P~ 
2 , 


U 


(8.71) 


v = 0 v “ " ' 

Substituting the value of p from Eq. (8.67) and then invoking the value of Q in 
Eq. (8.71), the final expression for shear stress becomes 


T 0 = 4 

n 


6flUh x h 2 


h\h i +h 2 ) 

The drag force required to move the lower surface at speed U is expressed by 


i 

D = J T 0 dx - 
o 


pUl 


h\ tb 


., hi hi h~) 

4In— - 6- 1 - 


fh 


\+h 1 


(8.72) 


Michell thrust bearing, named after A.G.M. Michell, works on the principles 
based on the theory of hydrodynamic lubrication. The journal bearing(Fig. 8.11) 
develops its force by the same action, except that the surfaces are curved. 



Fig. 8.11 Hydrodynamic action of a journal bearing 


8.5.2 Low Reynolds Number Flow Around a Sphere 


Stokes obtained the solution for the pressure and velocity field for the slow motion 
of a viscous fluid past a sphere. In his analysis, Stokes neglected the inertia terms 
of Navier-Stokes equations. Details of the solutions are beyond the scope of this 
text. However, integrating the pressure distribution and the shearing stress over 
the surface of a sphere of radius R, Stokes found that the drag D of the sphere, 
which is placed in a parallel stream of uniform velocity [/«,, is given by 

D = 6 Tip RU^ (8.73) 

This is the well-known Stokes’ equation for the drag of a sphere. It can be 
shown that one third of the total drag is due to pressure distribution and the 
remaining two third arises from frictional forces. If the drag coefficient is defined 
according to the relation 


C D 


D 

\ P UU 


(8.74) 
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where 



is the frontal area of the sphere, then 


or 


C D 

C D 


6n/uRU „ 


(8.75) 


A comparison between Stokes’ drag coefficient in Eq. (8.75) and experiments 
is shown in Fig. 8.12. The approximate solution due to Stokes’is valid for Re < 1. 



10~ 3 1(T 2 KT 1 1 10 10 2 10 3 10 4 10 5 10 6 Re 10 s 


Fig. 8.12 Comparison between Stokes' drag coefficient and experimental 
drag coefficient 


An important application of Stokes’ law is the determination of viscosity of a 
viscous fluid by measuring the terminal velocity of a fallingsphere. In this device, 
a sphere is dropped in a transparent cylinder containing the fluid under test. If the 
specific weight of the sphere is close to that of the liquid, the sphere will approach 
a small constant speed after being released in the fluid.Nowwecanapply Stokes’ 
law for steady creeping flow around a sphere where the drag force onthe sphere 
is given by Eq. (8.73). 

With the sphere, falling at a constant speed, the acceleration is zero. This 
signifies that the falling body has attained terminal velocity and we can say that 
the sum of the buoyant force and drag force is equal to weight of the body. 


y k r3 Psg = J nRi PiS + 6niiV T R 


(8.76) 


where p s is the density of the sphere, p, is density of the liquid and V T is the 
terminal velocity. 

Solving for fi, we get 


P = 


2 gR 2 

9 V T 


(Ps~ Pi) 


(8.77) 
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The terminal velocity V T can be measured by observing the time for the sphere to 
cross a known distance between two points after its acceleration has ceased. 

Summary 

• The Navier-Stokes equations, based on the conservation of momentum 
have been derived for a viscous incompressible fluid. 

• The Navier-Stokes equations are not amenable to an analytical solution 
due to the presence of nonlinear inertia terms in it. However, there are 
some special situations where the nonlinear inertia terms are reduced to 
zero. In such situations, exact solutions of the Navier-Stokes equations are 
obtainable. This includes the plane Poiseuille flow, the Couette flow, the 
flow through a straight pipe and the flow between two concentric rotating 
cylinders. All these flows are known as parallel flow where only one 
component of the velocity is non-trivial. 

• Knowledge of the velocity field obtained through analytical methods 
permits calculation of shear stress, pressure drop and flow rate. 
Applications of the parallel flow theory to the measurement of viscosity 
and hydrodynamics of bearing-lubrication are explained. 
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Solved Examples 


Example 8.1 _Two infinite plates are at h distance apart as in Fig. 8.13. There is a 

fluid of viscosity ft between the plates and the pressure is constant. The upper plate is 
moving at speed U — 4 m/s. The height of the channel h = 1.8 cm. Calculate the shear 
stress at the upper and lower walls if ft = 0.44 kg/m.s and p = 888 kg/m 3 . 



Fig. 8.13 Parallel flow between two plates with upper plate moving 

Solution Re = phU/fl = (888) (1.8/100) (4)/0.44 = 145. So, the flow is laminar and T = 

du . . , U 

ft—, u at any v is given by — y. 
oy ' h 

Shear stresses at the two walls are of equal magnitude, therefore, 
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T= =H (U ~ 0) (0.44)(4)/(l.8/100) 
dy h 

= 97.8 Pa 


Example 8.2 _ Water at 60° flows between two large flat plates. The lower plate 
moves to the left at a speed of 0.3 m/s. The plate spacing is 3 mm and the flow is laminar. 
Determine the pressure gradient required to produce zero net flow at a cross-section. (p = 
4.7 x 10” 4 Ns/m 2 at 60 °C) 


U = 0.3 m/sec 
b = 3 mm 


///////////////////////////// 


b 



Fig. 8.14 


_ , . _ d u dp 

Solution Governing equation: p — j = — 

dy dx 

1 dp 2 ^ „ 

u - — -— y + Qy + C 2 

2/u dr¬ 
at y = 0, u = -U, C 2 = -U 
at v = b, u = 0, which yields 

0= — ^-b 2 + Qb - U, 
111 dx 


or 


Cj = --—.*P-b 
b 2/u dv 


Now, 


u = ±.^ ( /- by)+ vU -1 


2 n dx 

b b p 


0= Jndy = J 


o 


0 


d|y_ w+t , [2_, 


dy 


or 


For, 


12 n dx 2 


Q = 0, with /i = 4.7 X 10 4 Ns/m' 


dp 6Up -6x0.3x4.7xl0” 4 n . nAT/ 2 
- - — =-94.0N/m-m 


dx 


(0.003) 


Example 8.3 A continuous belt (Fig. 8.15) passing upward through a chemical 
bath at velocity U 0 , picks up a liquid film of thickness, h, density p, and viscosity p. 
Gravity tends to make the liquid drain down, but the movement of the belt keep the fluid 
from running off completely. Assume that the flow is fully developed and that the 
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atmosphere produces no shear at the outer surface of the film. State clearly the boundary 
conditions to be satisfied by velocity at y = 0 and y = h. Obtain an expression for the 
velocity profile. 


Belt 



Solution The governing equation is 

A^u 

P~T^ = Pg 

d y 


or 


An , ^ 

p— = pgy + c \ 

Av 


d u = pgy + Q_ 

d y p p 

« = R f 1 + —y + c 2 

2 p p 

at y = 0, u - U 0 , so C 2 = U 0 
(\u 

at v = h, T = 0, so — =0 and C, = - pgh 

Ay 


= pgy~ _ pghy 

2 P P 


+ £/<>= — 


f v 2 ' 

1 - hy 

2 


+ U n 


Ex amp le 8.4 Water enters a rectangular duct at a rate of 10 m 3 /s as shown below. 
Two of the faces of the duct are porous. On the upper face, water is added at a rate shown 
by the parabolic curve, while on the front face water leaves at a rate determined linearly 
by the distance from the end. The maximum values of both rates are as given in Fig. 8.16. 
What is the average velocity leaving the duct if it is lm long and has a cross-section of 
0 .1m 2 ? 
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"2 

3 m /s per unit 



Solution Consider the top surface. The water enters the top surface in a parabolic 
manner. Let us first find this parabolic curve. 

Let w = ay 2 + by + c, where w is the flow rate per unit length at top face. Following are the 
boundary conditions: 
at v = 0, w = 0 
at y — 1, w = 3 


at 


,- 0.^-0 

d y 


So, c = b = 0, a = 3 
Thus, w = 3 y 2 

Similarly, consider the front surface, 
let u = my + d i 

when y= 1,« = 0 r we get: d = 5 

y = 0, u = 5 m — —5 

So, u = - 5y + 5 

For steady incompressible flow, the continuity equation gives 


J V-dA =0 

CS 

Choose the interior of the duct as a control volume. 


l l 

Thus, - 10- J* 3j^ 2 ciy + J(5-5y)dy+ ^(O.l) =0 

o o 

or -10-l + f5-|j =-0.1 V 2 
Finally, V 2 = 85 m/s 


Ex amp le 8.5 Water at 20 °C is flowing between a two dimensional channel in 
which the top and bottom walls are 1.5 mm apart. If the average velocity is 2 m/s, find out 
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(a) the maximum velocity, (b) the pressure drop, (c) the wall shearing stress, and (d) the 
friction coefficient [p = 0.00101 kg/m ■ s]. 

Solution (a) The maximum velocity is given by 
t / max= | t/ av = |[ 2 ] = 3 m/s 

(b) The pressure drop in a two dimensional straight channel is given by 

— = — ^ ^ lnax , where b = half the channel height 

dx b~ 

-2(0.00101)3 

[1.5/(2 x 1000)] 2 

= - 10773.33 N/m 3 

or -10773.33 N/m 2 per metre 

(c) The wall shearing stress in a channel flow is given by 


du , d» —1.5 

x - b — =- 

y dy dx- 2x1000 

or x yx — 8.080 N/m 2 

(d) Re (based on channel height and average velocity) 


(-10773.33) 


pU m (2b) 1000 x 2 x (1,5/1000) otV7A 
p 0.00101 


which is more than 2300 but the flow is not turbulent as well. Laminar approximation is 


quite logical. However, C7= 


12 

Re 


0.004. 


Example 8.6 The analysis of a fully developed laminar flow through a pipe can 
alternatively be derived from control volume approach. Derive the expression 


v- = 


R 

4/7 


dz 


"7 


starting from the control volume approach. 


dr 




Fig. 8.17 Fully developed laminar flow through a pipe 


Solution Let us have a look at Fig. 8.17. The fluid moves due to the pressure gradient 
which acts in the direction of the axis and in the sections perpendicular to it the pressure 
may be regarded as constant. Due to viscous friction, individual layers act on each other 


producing a shearing stress which is proportional to 


dv-_ 

dr 





























The McGraw-Hill Companies 



Viscous Incompressible Flows 


1^47 


In order to establish the condition of equilibrium, we consider a fluid cylinder of length 
51 and radius r. Now we can write 

[p - (p + 8p)] k r = - t 2 n r 51 

or - 8p Jir 2 = - id 2 nr 81 

dr 


upon integration, 


dv, _ 1 dp _ 1 dp 

dr Ijj. dl 2p dr 


1 dp i 

=-— + K 


4p dr 


at r =R, v „ = 0, hence K = - 


J_ dp 
4p dz 


■R* 


So, = 

4/u V dr 




In the laminar flow of a fluid in a circular pipe, the velocity profile is 
The rate of discharge is then represented by volume of a paraboloid. 
Prove that for this case the ratio of the maximum velocity to mean velocity is 2. 
Solution See Fig. 8.17. For a paraboloid, 



Example 8.7 

exactly a parabola. 


Q = 


R 

( r\ 2 

& 

II 

1- - 

J J max 

UJ 

0 



(2k r dr) 


= 2k v. 


2 4 R 2 


= 2 KV- 


= 


( kR 2 ^ 


V 2 , 


Q V- (kR" 12.) V. 

U = _ max _ _ max 

Zmean A~ (kR 2 ) ~ 2 



Thus, 


= 2 


Example 8.8 _ The velocity distribution for a fully-developed laminar flow in a pipe 
is given by 
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u = ^.?P [l - {r/R f ] 

4p dz 

Determine the radial distance from the pipe axis at which the velocity equals the 
average velocity. 

Solution For a fully-developed laminar flow in a pipe, we can write 


Rr_ dp_ 
4p dz 


-'i 


V = — 

av A n 

A nR 


_ i 

kR 2 J 1 4 p dz 


= R^dp 

&p dz 

Now, for u =K we have, 


B?_ dp_ 
4p dz 


1- 


B?_ dp_ 
'ip dz 


'-'S 


• 2nrdr 


or 


-i 


r = -£=• = 0.707 R 

V2 


Exam ple 8.9 _ SAE 10 oil is flowing through a pipe line at a velocity of 1.0 m/s. The 
pipe is 45 m long and has a diameter of 150 mm. Find the head loss due to friction, [p = 
869 kg/m 3 and p = 0.0814 kg/m ■ s] 

yi 


Solution 


h f = 


fiv 1 

2gD 


In order to know f first we have to calculate Re. 
Re= 


(869) (1) (150/1000) 


= 1601.35 


0.0814 

Since Re < 2000, the flow can be assumed to be laminar and 
/= 64/Re = 64/1601 =0.04 


So, 


h r 


(0.04)(45)(1.0) 2 
(2)(9.81)(150/1000) 


= 0.612 m 
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Example 8.10 Heavy fuel oil flows from A to B through a 100 m horizontal steel 
pipe of 150 mm diameter. The pressure at A is 1.08 MPa and at B is 0.95 MPa. The 
kinematic viscosity is 412.5 x 10 -6 m 2 /s, and the relative density of the oil is 0.918. What 
is the flow rate in m 3 /s? 


Solution The Bernoulli’s equation between A and B 


or 


1.08 xlO 6 V 2 ,, 

-1-h 0 

918x9.81 2 g 


0.95 xlO 6 
918x9.81 


V 

+ — + 0 + / 
2g 


100F 2 
2g x 0.150 


119.925-105.49 = 


666.67 fV 2 
2 g 


Both V and /are unknown so we have to follow another approach. If laminar flow exists, 
then from Eq. (8.48), 


or 


Pi-Pi = 
l R 2 

v = (Pi-Pi)D 2 _ (1080-950) 1000(.150) 2 
av 32 fil ~ 32(412.5xl0~ 6 x918)(100) 


130x1000x0.0225 

32x412.5x918x100 


= 2.41 m/s 


Re = 2,41 x (150/1000) =g7636 
412.5x10 


Hence, laminar flow assumption is valid. 


Q = AV W = -^ (0.15) 2 x 2.41 =0.0425 m 3 /s 


Example 8.11 A wind tunnel has a wooden (e = 0.0001 m) rectangular section 40 
cm by 1 m by 50 m long. The average velocity is 45 m/s for air at sea-level standard 
conditions. Find the power required if the fan has 65 percent efficiency. For air, p = 
1.2 kg/m 3 , p = 1.81 x 10~ 5 kg/m-s. 

Solution The wetted perimeter of the duct is 


= 


40 40 , , 

-1-h 1 + 1 

100 100 


= 2.8 m 


40 2 

and the flow area is =-x 1 = 0.4 m . 

w 100 


Hence, hydraulic diameter D h = 


4x0.4 

2.8 


= 0.5714 m 


Re = 


e/D h = 


P Dh Ktv 

P 

0.0001 

0.5714 


_ (1.20)(0.5714) (45) 

(1.81 x 10 5 ) 

0.000175 


= 1.7x 10 6 


Now, 
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From Moody’s chart for Re = 1.7 x 10 6 and eJD h = 0.000175 
/= 0.0140 


h f=f 


V Ai y 


V 

— = (0.0140) 
2 g 


f 50 1 

" (45) 2 " 

[ 0.5714 J 

2x9.8 


= 126.5 m 

and pressure drop, A p = pgh f = (1.20) (9.8) (126.5) = 1489 Pa 


Pumping power, 


e= ^ F= (S (i)(45)=i8 ' 0m3/s 

PgQhf _ (1.20)(9.8)(18.0) (126.5) 


n 

= 41200W 


0.65 


Example 8.12 A circular pipe of radius a and length L is attached to a smoothly 
rounded outlet of a liquid reservoir by means of flanges and bolts as shown in 
Fig. 8.18. At the flange section the velocity is uniform over the cross-section with 
magnitude V 0 . At the outlet, which discharges into the atmosphere, the velocity profile is 
parabolic because of the friction in the pipe. What force must be supplied by the bolts to 
hold the pipe in place? 


F CV 



Solution For the control volume as shown, continuity equation gives (for steady flow) 

a a 2 

jpVdA=0 or, - V 0 ncc + J2^;rrdr = 0, So, J V e rdr= ^° a 
cs oo 

since V e is parabolic, let V e — a 0 + ap- + a 2 r 2 
at r = a,V e = 0 —» a x a + a 2 a 2 + a 0 = 0 
dV 

at r = 0, —i = 0 a, = 0 
dr 1 

This will give 



2x , V 0 a 1 

+ a 2 r )rdr = —— 


o 


2 
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or 


2 4 2 

a a Tr a 

a n -b a 2 — = V n — 

u 2 1 4 u 2 

2 


fl 0 + ? a 2 Vo 


Again, we know a 0 + a a 2 = 0 

Combining the above two expressions, we get, 

2 


— a 2 = - V 0 a 2 = - '—t~ and a 0 = + 2 V 0 
2 a 


V= 2Vr 


2 ' 

1-^T 

a 


Now, applying momentum equation to the control volume in the flow direction (let F be 
the force as shown on the control volume; same force must be supplied by bolts): 


u 

F + (Pi -Pi) Kal = - PV o na 1 ■ V 0 + p J 2jirdr V} 


= - up Vo a 2 + 


-KpVo J ! 


.2 A 

r "-v 


dr 


2 2 


= K p V 0 a 


= np Vo a 2 


1 

, f r 

( 4 

r „ r 

f r ^ 

-1+ 8- 

1 + —4 ~ 2 ~2 

d - 

J a 

a q a z 

\a J 

0 

\ / 



-1 + 1 4 + --4 


= ^npa 2 V 0 2 

1 2 2 2 

F = — npar V 0 - (p t - p-^na in horizontal direction only (gravity is in vertical 
direction). 

Example 8.13 A slipper (slider) and plate (guide), both 0.5 m wide constitutes 
a bearing as shown in Fig. 8.19. Density of the fluid, p = 9.00 kg/m J and viscosity, p = 
0.1 Ns/m 2 . 



Fig. 8.19 Slipper and plate both 0.5 m wide 

Find out the (a) load carrying capacity of the bearing, (b) drag, and (c) power lost in 
the bearing. 
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Solution (a) Considering the width as b and using Eq. (8.70) for the load carrying 
capacity, 


P = 


6nUl A b 


ln k_ 2 hzh. 

/? 2 A, + h 2 


Oh ~ ^2 ) 

6 x 0.1 x 1 x 0.2 x 0.2 x 0.5 
(0.005 - 0.002) 2 

" 0.005 0.005 -0.002 

11 0.002 _ 0.005 + 0.002 

0.012 


^ (0.9163-0.8571) = 78.93 N 
9.0 xlO -6 

(b) Making use of Eq. (8.72) for width b, the drag force may be written as 

h ~ h 2 


D= HUlb 


\ — h 2 


4 In— - 
lb 


0.1 x 1 x 0.2 x 0.5 


(0.005-0.002) 
0.01 


\ + !h 


0.005 0.005 - 0.002 

4 In-6 


0.003 


0.002 0.005 + 0.002 

(3.6651 -2.5714) = 3.645 N 


(c) Power lost = drag x velocity 

= 3.645 x 1 =3.645 W 


Exam pie 8.14 A cylindri¬ 
cal journal bearing supports a 
load directed vertically up¬ 
wards, with the shaft rotating 
clockwise. Sketch the position 
of the shaft centre with respect 
to that of the bearing (hole), if 
no cavitation is present. Give 
explanation. No equations are 
required. 

Solution In this case, the 
pressure distribution is sym¬ 
metric as shown (Fig 8.20) 
where 6 is measured in the 
direction of the rotation of shaft, 
from the position of maximum 
clearance. Thus, as shown, the 
shaft centre is to the right of 
bearing centre, and the line of 
centres is horizontal. 



(a) 
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Example 8.15 For the following thrust bearing (Fig. 8.21), show that the force on 
the straight slider in the x-direction is the same as that on the guide. 


Po 



- 


p 0 = atm pressure 

Fig. 8.21 

Solution It is given that the velocity profile is 


and load 


u 

U 


P = 


fl-^1 

1 - 3 — 

( 

1- 

2 

h Y 

l hj 

h 

V 

n + 1 

h)\ 


6pUL z 




In n ■ 


2(»-l) 

n +1 


where n = h l /h 2 

Force on the slider in the x-direction is 

L 


F = 


f , cos a f, . dx ... . 

T s dx(l)-+ (/> “ Po) -(1) sin 

J cos a j cos a 

o o 

L L 

Jtj dx + tan a j ( ^ -p 0 )dx 


a 


Now, 


du 


' v — h 


and 

m= t/fl--l 

1 — 3 — 

( 

1- 

2 

h V 


l h) 

h 

Y 

n + 1 

h) 


, so we get 


du 

dv 


— =U 


—— — 3 I — — —y 
h h 


,_2_A 

» + l h 


T s =-pU 


1 3 

-1- 

h h 


1- 


2 h 


n + I h 
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Also, 


Thus, 


= llU 


2 6 h 

-1-.—y 

h n +1 h~ 


d h = 


T. dx = 


/ ^dx 

L 

T h C 

J d h 


h, - h 


h = h ] - (h x - h 2 ) — 


*1 , 

Atzc/ fr 

2 6 /z. 

h -h 2 J l 

h n + 1 h 2 , 


Ah 


Also, load 


pUL 


h\ - h 2 
HUL 


2 In h ■ 


6 /?] 


h 2 (n — 1) 
\lUL 


P = 


h 2 (n-l) 

L 

J (P~Po) 


2 In n — 


- 2 In n + 


n +1 h 
6 1\ 


n +1 

6 (n -1) 
n +1 



--T 

A 

*2 


dx cos a 


cos a 


(neglecting contribution of 
T s to load; a is small) 


F s = J r s dx + tan a ( P ) 


F = 


\iUL 


h 2 {n -1) 
6/J.UL 2 


2 In n + 


6(n-\) 
n +1 



Ai - 

J 

Z, 


h 2 (n -1 ) 2 
HUL 


h 2 (n -1) 

Now the force on the guide is 

L 

F r . = \r n dx 


In n ■ 


4 In n- 


2(»-l) 
« +1 

6(n -1) 
n +1 


G~ J ‘-G'- 
o 

%= -F 


du ' 

^vj v: 


/y=0 


But 
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= n u 


The expression is T G = /J.U 


1 3 

-1- 

h h 


\ _2_A' 

n + 1 h 


6 /?! 


h n + 1 h 1 


F g = r G dx 


/jj — h 2 


"i 

Jigd/; (as before) 


\lUL 4 6 /j 




h 2 (n — \)j\h 77 + 1 h y 


d/? 


[lUL 


^2 (» - 1 ) 
luUL 


A i 6/7] 

4 In 77 +-L 


77 + 1 


/? 2 ( 77 - 1 ) 


4 In 77 - - 


77+1 


( 1 

i y 

A 

h 2 ) 

d" 



which is the same as /+. 


Exercises 


8.1 Choose the correct answer. 

(i) Bulk stress is equal to thermodynamic pressure 

(a) if second coefficient of viscosity is zero 

(b) for incompressible flows 

(c) for a compressible fluid with negligible second coefficient of viscosity 

(d) if bulk coefficient of viscosity is non-zero. 

(ii) Assumptions made in derivation of Navier-Stokes equations are: 

(a) continuum, incompressible flow Newtonian fluid and fJ. = constant 

(b) steady flow, incompressible flow, irrotational flow 

(c) continuum, non-Newtonian fluid, incompressible flow 

(d) continuum, Newtonian fluid, Stokes’ hypothesis and isotropy. 

(iii) In a fully developed pipe flow 

(a) pressure gradient is greater than the wall shear stress 

(b) inertia force balances the wall shear stress 

(c) pressure gradient balances the wall shear stress only and has a constant 
value 

(d) none of the above 

(iv) In the case of fully developed flow through tubes, 

(a) Darcy’s friction factor is four times the skin friction coefficient 

(b) Darcy’s friction factor and skin friction coefficients are same 

(c) Darcy’s friction factor is double the skin friction coefficient 
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8.2 


(d) the skin friction coefficient is greater than the Darcy’s friction factor, 
(v) Based on hydrodynamic theory of lubrication, state which of the following 
are correct. 

(a) The load bearing capacity remains unchanged so long either the slipper 
or the bearing moves in the same direction while the other is held fixed. 

(b) Reversing the direction of the movement of the slipper, bearing 
remaining fixed, does not cause any change in load bearing capacity. 


, . du d 2 u 

(c) u—»n— T 

ox oy 

(d) For a large film thickness, h{x), the maximum pressure location shifts 
from the middle. 

(vi) Observation on a spherical object falling in a liquid pool is the method of 
measuring viscosity by making use of Stokes’ viscosity law. The falling 
body attains terminal velocity if 

(a) the weight of the falling body is more than the sum of the buoyancy 
force and the drag force 

(b) the drag force is equal to the buoyancy force 

(c) the buoyancy force is more than the drag force 

(d) the sum of the buoyancy force and the drag force is equal to the weight 
of the body. 

(a) What is the basic difference between the Euler’s equations of motion and the 
Navier-Stokes equations? 

(b) In case of flow through a straight tube of circular cross-section with 
rotational symmetry, the axial component of velocity is the only non-trivial 
component and all the fluid particles move in the same direction only. Find 
out the average velocity and the maximum velocity within the tube. If Darcy- 
Weisbach equation for pressure drop over a finite length is given by h,=f (LI 
D ) (V 2 /2g), prove that/= 64/Re, where L is the length andD is the diameter 
of the tube. 


8.3 What is the relationship between the average velocity and maximum velocity in 
case of parallel flow between two fixed parallel plates? What do you understand 
by inlet region and developed region? 

Ans. (C/ max =1.5 C/ av ) 

8.4 Show that in case of a Couette flow, the shear stress at the horizontal mid-plane 
of the channel is independent of the pressure gradient imposed on the flow. 

8.5 (a) Find out the total load and the frictional resistance on a block moving with a 

velocity U over a horizontal plate separated by a thin layer of lubricating oil, 
the thickness of layer being h l and h 2 at the edges of the block which has a 
straight bottom. 

(b) Also show that the volume flow rate of lubricant is given by 


< 2 = U 


hi + h-, 


8.6 Oil flows between two parallel plates, one of which is at rest and the other moves 
with a velocity U, (a) If the pressure is decreasing in the direction of flow at the 
rate of 5 Pa/m, the dynamic viscosity is 0.05 kg/ms, the spacing of the horizontal 
plate is 0.04 m and the volumetric flow Q per unit width is 0.02 m 2 /s, what is the 
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8.7 


8.9 


8.10 


8.11 


8.12 


8.13 


velocity W. (b) Calculate U if the pressure is increasing at a rate of 5 Pa/m in the 
direction of flow. 

Ans. (a) 0.97 m/s (b) = 1.027 m/s 
Water flows between two very large, horizontal, parallel flat plates 20 mm apart. 
If the average velocity of water is 0.15 m/s, what is the shear stress (a) at the 
lower plate, and (b) 5 mm and 10 mm above the lower plate? Assume p = 1.1 x 


10~ 3 Ns/m 2 . 


Ans. 


(a) 0.0495 N/m 2 (b) 0.0248 N/m 2 


A Newtonian liquid flows slowly under gravity along an inclined flat surface that 
makes an angle 6 with the horizontal plane. The film thickness is T and it is 
constant. The flow is two-dimensional, (a) Show that the fluid velocity u along x 
(flow) direction is given by 


g sin 6 


T- — 


(b) 


v V 2, 

Calculate the average velocity w av , and the volumetric flow rate Q per unit 
width of the surface. The pressure within the fluid is a function of y alone, 
where y is the nonnal to the flow direction. The ^-component of velocity is 
trivial. 

Ans. (I/ av = g sin 6 T l l'iv, Q = g sin 

A horizontal circular pipe of outer radius R { is placed concentrically inside 
another circular pipe of inner radius R 2 . Considering fully developed laminar 
flow in the annular space between pipes show that the maximum velocity occurs 
at a radius R 0 given by 

-il/2 


Ro = 


2 In (R 2 IR { ) 


The Reynolds number for flow of oil through a 5 cm diameter pipe is 1700. 
The kinematic viscosity, v= 1.02 x 10~ 6 m 2 /s. What is the velocity at a point 
0.625 cm away from the wall. 

Ans. (0.03 m/s) 

The velocity along the centre line of the Hagen-Poiseuille flow in a 0.1 m diameter 
pipe is 2 m/s. If the viscosity of the fluid is 0.07 kg/ms and its specific gravity is 
0.92, calculate (a) the volumetric flow rate, (b) shear stress of the fluid at the pipe 
wall, (c) local skin friction coefficient, and (d) the Darcy friction coefficient. 

Ans. (a) 7.854 x 10~ 3 m 3 /s (b) 5.6 N/m 2 , (c) 0.012 (d) 0.048 
Kerosene at 10 °C flows steadily at 20 1/min through a 150 m long horizontal 
length of 5.5 cm diameter cast iron pipe. Compare the pressure drop of the 
kerosene flow with that of the same flow rate of benzene at 10 °C through the 
same pipe. For kerosene at 10 °C, p = 820 kg/m 3 and p = 0.0025 Ns/m 2 and for 
benzene p = 899 kg/m 3 and p = 0.0008 Ns/m 2 . Why do you obtain greater 
pressure drop for benzene? 

A viscous oil flows steadily between parallel plates. The fully developed velocity 
profile is given by 


h 2 


i 

(Ml 

8p 



l h ) 
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where the total gap between the plates is h = 3 mm andy is the distance from the 
centre line. The viscosity of the oil is 0.5 Ns/m 2 and the pressure gradient is 
- 1200 N/nr/m. Find the magnitude and direction of the shear stress on the upper 
plate, and the volumetric flow rate per metre width of the channel. 

Ans. (a) -1.80 N/m 2 , (b) 5.40 x 10~ 6 m 3 /s m 
8.14 A fully developed laminar flow is taking place in the annulus between two 
concentric pipes. The inner pipe is stationary, and the outer pipe is moving in the 
axial direction with a velocity V 0 . Assume the axial pressure gradient to be zero 
(dy>/dz = 0). Find out a general expression for the shear stress as a function of 
radial coordinate. Also find out a general expression for the velocity profile VJr). 

Ans. (a) T = Air, (b) V z = V 0 


In (r 0 /r,) 
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Laminar Boundary Layers 


9.1 INTRODUCTION 

The boundary layer of a flowing fluid is the thin layer close to the wall. In a flow 
field, viscous stresses are very prominent within this layer. Although the layer is 
thin, it is very important to know the details of flow within it. The main-flow 
velocity within this layer tends to zero while approaching the wall. Also the 
gradient of this velocity component in a direction normal to the surface is large as 
compared to the gradient of this component in the streamwise direction. 

9.2 BOUNDARY LAYER EQUATIONS 

In 1904, Ludwig Prandtl, the well known German scientist, introduced the 
concept of boundary layer [1] and derived the equations for boundary layer flow 
by correct reduction of Navier-Stokes equations. He hypothesized that for fluids 
having relatively small viscosity, the effect of internal friction in the fluid is 
significant only in a narrow region surrounding solid boundaries or bodies over 
which the fluid flows. Thus, close to the body is the boundary layer where shear 
stresses exert an increasingly larger effect on the fluid as one moves from free 
stream towards the solid boundary. However, outside the boundary layer where 
the effect of the shear stresses on the flow is small compared to values inside the 
boundary layer (since the velocity gradient du/dy is negligible), the fluid particles 
experience no vorticity, and therefore, the flow is similar to a potential flow. 
Hence, the surface at the boundary layer interface is a rather fictitious one 
dividing rotational and irrotational flow. PrandtTs model regarding the boundary 
layer flow is shown in Fig. 9.1. Hence with the exception of the immediate vicinity 
of the surface, the flow is frictionless (inviscid) and the velocity is U. In the 
region, very near to the surface (in the thin layer), there is friction in the flow 
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which signifies that the fluid is retarded until it adheres to the surface. The 
transition of the mainstream velocity from zero at the surface to full magnitude 
takes place across the boundary layer. Its thickness is <5 which is a function of the 
coordinate direction x. The thickness is considered to be very small compared to 
the characteristic length L of the domain. In the normal direction, within the thin 
layer, the gradient du/dy is very large compared to the gradient in the flow 
direction duldx. Next step is to simplify the Navier-Stokes equations for steady 
two dimensional laminar incompressible flows. Considering the Navier-Stokes 
equations together with the equation of continuity, the following dimensional 
form is obtained. 


du du 

U lh +V lty 


dv dv 

u - v v — 

dx dy 


1 dp 
p dx 

1 dp 

P dy 


P 

P 

a 

p 


d 2 u 

dx 2 

d 2 v 

dx 2 


d 2 u 

df 

d 2 v 

dy 2 


(9.1) 

(9.2) 


du + dv _ o 
dx dy 


(9.3) 


t/« 



Fig. 9.1 Boundary layer on a flat plate 


Here the velocity components u and v are acting along the streamwise x and 
normal y directions respectively. The static pressure is p, while p is the density 
and ,U is the dynamic viscosity of the fluid. 

The equations are now non-dimensionalised. The length and the velocity 
scales are chosen as L and U m respectively. The non-dimensional variables are: 


* 

u 




p 


pul 


where [/«, is the dimensional free stream velocity and the pressure is non- 
dimensionalised by twice the dynamic pressure p d = (1/2) p U 2 m . Using these 
nondimensional variables, the Eqs (9.1) to (9.3) become 







































The McGraw-Hill Companies 




Laminar Boundary Layers 


OH 


* 

u 




* 

u 






1 

Re 


d 2 u 
dx 2 


1 

Re 


d 2 v 

dx* 2 


n2 * 

a u 

-i2 * 

a v 


(9.4) 

(9.5) 


du dv 


= 0 


dx dy 
where the Reynolds number, 

pU x 


(9.6) 


L 


Re = 

P 

Let us examine what happens to the u velocity as we go across the boundary 
layer. At the wall the u velocity is zero. The value of u on the inviscid side, that is 
on the free stream side beyond the boundary layer is U. For the case of external 
flow over a flat plate, this U is equal to U„. 

Based on the above, we can identify the following scales for the boundary 
layer variables: 


Variable 

Dimensional scale 

Nondimensional scale 

u 

(A 

1 

X 

L 

1 

y 

5 

e(=S/L) 


The symbol e describes a value much smaller than 1. Now, let us look at the 
order of magnitude of each individual terms involved inEqs(9.4), (9.5) and (9.6). 
We start with the continuity Eq. (9.6). One general rule of incompressible fluid 
mechanics is that we are not allowed to drop any term from the continuity 
equation. From the scales of boundary layer variables, the derivative du /dx is of 
the order 1. The second term in the continuity equation dv /dy should also be of 

_ * * _ x- 

the order 1. Now, what makes dv /dy to have the order 1? Admittedly v has to 
be of the order e because y becomes e (= S/L) at its maximum. Next, consider 
Eq. (9.4). Inertia terms are of the order 1. Among the second order derivatives, 
d 2 u/dx* 2 is of the order 1 and d 2 u*/dy* 2 contains a large estimate of (1/e 2 ). 
However after multiplication with 1/Re, the sum of these two second order 
derivatives should produce at least one term which is of the same order of 
magnitude as the inertia terms. This is possible only if the Reynolds number (Re) 
is of the order of 1/e 2 . It follows from the Eq. (9.4) that -dp Idx will not exceed 
the order of 1 so as to be in balance with the remaining terms. Finally, Eqs (9.4), 
(9.5) and (9.6) can be rewritten as 


* du* 

* du* 

V dy 

dp* 

dx* 

1 

2)2 * 2)2 * 

d u d u 

u * 
dx 

Re 

-> *2 -> *2 
dx dy 

(1) 

<«>£ - 
(e) 

(1) 

< 4 ®+ ' 1 
(1) (£ 2 ) 


(9.4) 
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dv* 


dv* 


dp 


+ 


d\ 


d\ 


dx* 

dy* 

dy* 

Re 

dx* 2 

a )tn 
0 ) 

(£) 77 ‘ 
(£) 

(?) 

(£ 2 ) 

(e) 

(!) 

du* 

dv* 

dy* 

0 



dx* 



(1) 

(£) 




(1) 

(£) 





*2 


dy 


(£ 2 ) 


(9.5) 


As a consequence of the order of magnitude analysis, 


dx 


*2 


(9.6) 


can be dropped 


from the x direction momentum equation, because on multiplication with 1/Re it 
assumes the smallest order of magnitude. Now, consider the y direction momen¬ 
tum Eq. (9.5). All the terms of this equation are of a smaller magnitude than those 
of Eq. (9.4). This equation can only be balanced if dp !dy is of the same order of 
magnitude as other terms. Thus they momentum equation reduces to 


14 = O (£) 

dy 


(9.7) 


This means that the pressure across the boundary layer does not change. The 
pressure is impressed on the boundary layer, and its value is determined by 
hydrodynamic considerations. This also implies that the pressure p is only a 
function of x. The pressure forces on a body are solely determined by the 
inviscid flow outside the boundary layer. The application of Eq. (9.4) at the 
outer edge of boundary layer gives 

* d u* dp* 

U -— =-=r 

dx dx 


(9.8a) 


In dimensional form, this can be written as 

dU = _1_ dp 

dr p dr 

On integrating Eq. (9.8b), the well known Bernoulli’s equation is obtained, 


U 


(9.8b) 


1 2 

— p U = a constant 


(9.9) 


Finally, it can be said that by the order of magnitude analysis, the Navier-Stolces 
equations are simplified into equations given below. 


du 

dx* 


du* _ dp* 
dy* " 

dp* 


1 d 2 


dx 


Re dy 


*2 


dli* 

dx* 


dy* 

dv* 

dy* 


= 0 


= 0 


(9.10) 


(9.11) 


(9.12) 
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These are known as Prandtl’s boundary-layer equations. The available 
boundary conditions are: 

Solid surface 

aty = 0, u = 0 = v 
or at y = 0 u = 0 = v 

Outer edge of boundary-layer 

at y* = (e) = , u = 1 

or at y = 8, u = U (x) 

The unknown pressure p in the x-momentum equation can be determined 
from Bernoulli’s Eq. (9.9), if the inviscid velocity distribution U (x) is also 
known. The preceding derivations are related to a flat surface, but these can be 
easily extended to curved surfaces. While doing so, it is seen that Eqs (9.10) to 
(9.14) continue to be applicable only if the curvature does not change abruptly. 
Elowever, the boundary layer equations are relatively easier to solve as com¬ 
pared to the Navier-Stolces equations and have been solved by various analytical 
and numerical techniques. 

We solve the Prandtl boundary layer equations for u (x, y) and v (x, y) with 
U obtained from the outer inviscid flow analysis. The equations are solved by 
commencing at the leading edge of the body and moving downstream to the 
desired location. Note that the reduced momentum Eq. (9.10) is still nonlinear. 
However, it does allow the no-slip boundary condition to be satisfied which 
constitutes a significant improvement over the potential flow analysis while 
solving real fluid flow problems. The Prandtl boundary layer equations are thus 
a simplification of the Navier-Stokes equations. 

9.3 BLASIUS FLOW OVER A FLAT PLATE 


(9.13) 

(9.14) 


The classical problem considered by H. Blasius was a two-dimensional, steady, 
incompressible flow over a flat plate at zero angle of incidence with respect to 
the uniform stream of velocity U m . The fluid extends to infinity in all directions 
from the plate. The physical problem is already illustrated in Fig. 9.1. 

Blasius wanted to determine (a) the velocity field solely within the boundary 
layer, (b) the boundary layer thickness (8), (c) the shear stress distribution on 
the plate, and (d) the drag force on the plate. 

The Prandtl boundary layer equations in the case under consideration are 


du du 


dx 


dy 


= v 


du + dv _ o 
dx dy 


d 2 u 

df 

(9.3) 


The boundary conditions are 
aty = 0, 
at y = 00 


(9.15) 


u = v = 0 
u=U 


(9.16) 
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It may be mentioned that the substitution of the term 


J_ d p 

p dx 


in the original 


boundary layer momentum equation in terms of the free stream velocity produces 
d U 1 

U^ —— which is equal to zero. Hence the governing Eq. (9.15) does not 
dx 


contain any pressure-gradient term. However, the characteristic parameters of 
this problem are U M V, x, y, that is, 

u = u (U M V, x, y) 


Before we write down this relationship in terms of two non-dimensional 
parameters, we have to be acquainted with the law of similarity in boundary layer 
flows. It states that the u component of velocity with two velocity profiles of 
u(x, y) at different x locations differ only by scale factors in u and y. Therefore, 
the velocity profiles m(x, y) at all values of x can be made congruent if they are 
plotted in coordinates which have been made dimensionless with reference to the 
scale factors. The local free stream velocity U(x) at section x is an obvious scale 
factor for u, because the dimensionless u(x) varies between zero and unity with y 
at all sections. The scale factor fory, denoted byg(x), is proportional to the local 
boundary layer thickness so that y itself varies between zero and unity. The 
principle of similarity demands that the velocity at two arbitrary x locations, 
namely, x 1 and x 2 should satisfy the equation 


u\xx,{ylg{x i)}] = »[x 2 ,{y/g(x 2 )}] 

U(x i) U(x 2 ) 

Now, for Blasius flow, it is possible to write 

/ \ 


(9.17) 


U„ 


= F 


y 




= F{r]) 


(9.18) 


where p~ — and 5~ jdFF 


vx 


or more precisely, t] = 


y 


vx 

u 


(9.19) 


The stream function can now be obtained in terms of the velocity components 


as 


y/= J»dy = jV„ F(ri) dp = vx J F(rj) dr? 


or y/= vx ftp) + constant (9.20) 

where J F(p) dp = ftp) and the constant of integration is zero if the stream 
function at the solid surface is set equal to zero. 
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Now, the velocity components and their derivatives are: 

L=aX-.*L = u _ r{tf) 

dy drj dy 


v = 


'dy + dy dr]' 

v dx d r] dx y 




2 \!x 


1 


7 


(9.21a) 


1 


2 jvx x 

lu. 


or 


* = \r\f\m-fh 7)] 


(9.21b) 


or 


or 


du 

dx 

du 

dx 

du 

dy 

du 

dy 


u-f"0i)^-=u„f"(T 1 y 

dx 


—r~‘— 
2 x 


1 


y 


l 


2 jvx x 

i[7 


(9.21c) 


dy 


1 


= U„ 


u 


vx 


f"(v) 


(9.21d) 


d 2 


or 


dy 


d' i 

~d? 


U = 

vx 


1 


ui 


VX 


f" 0?) 


(9.21e) 


Substituting (9.21) into (9.15), we have 

7-2 


ui 


u ~ f'(j]) /"in) + ^=4n/'(n) - f(n)\ f"(n ) = ^ / w, (n) 


2 x 


2x 


or 


i ui 


-ir—mrxvi). 

2 x 


Ui 


rxn) 


or 2f" ’ (r\) + j{r]) f”(r\) = 0 (9.22) 

This is known as Blasius Equation. The boundary conditions as in Eq. (9.16), in 
combination with Eq. (9.21a) and (9.21b) become 
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at 77 = 0 : 7) = 0 , fin) = 0 

at n = °o : fin) = 1 

Equation (9.22) is a third order nonlinear differential equation. Blasius 
obtained the solution of this equation in the form of series expansion through 
analytical techniques which is beyond the scope of this text. However, we shall 
discuss a numerical technique to solve the aforesaid equation which can be 
understood rather easily. 

It is to be observed that the equation for/does not contain x. Further boundary 
conditions at x = 0 andy = °° merge into the condition q —> 00 , ulU„ =/' = 1. This 
is the key feature of similarity solution. 

We can rewrite Eq. (9.22) as three first order differential equations in the 
following way 

f' = G (9.24a) 

G' = H (9.24b) 

H'=-±fH (9.24c) 

Let us next consider the boundary conditions. The condition /(0) = 0 remains 
valid. Next the condition/'(0) = 0 means that G(0) = 0. Finally/'(°°) = 1 gives us 
G(oo) = l. Note that the equations for/and G have initial values. However, the 
value for//(0) is not known. Hence, we do not have a usual initial-value problem. 
Nevertheless, we handle this problem as an initial-value problem by choosing 
values of Hi 0) and solving by numerical methods ft 7 ), Gin), and H(n). In 
general, the condition G(°°) = 1 will not be satisfied for the function G arising 
from the numerical solution. We then choose other initial values of H so that 
eventually we find an //(0) which results in G(°°) = 1. This method is called the 
shooting technique. 

In Eq. (9.24), the primes refer to differentiation wrt the similarity variable 77 . 
The integration steps following Runge-Kutta method are given below 

fn+i =fn + ' (*. + 2 k 2 + 2k 3 + k 4 ) (9.25a) 

o 

G n+ 1 = G„+ ' (/, + 2/ 2 + 24 + / 4 ) (9.25b) 

o 

H n+l = H n + — (w 7 j + 2 m 2 + 2 m 3 + m 4 ) (9.25c) 

6 

One moves from q n to V n +i = n n + h- A fourth order accuracy is preserved if h 
is constant along the integration path, that is, q n +\ ~ r ln = h for all values of n. The 
values of k, l and m are as follows. 

For generality let the system of governing equations be 

f = F x if, G, H, 77), G' = F 2 if G, H, 77) and H' = F 3 if G, H, 77). 


(9.23) 


Then, 


K=hF x if G n , H n , n„) 
1\ = h F 2 if G n , H n , n„) 
— h F 2 if n , G n , H n , 1)n) 






The McGraw-Hill Companies 




Laminar Boundary Layers 


067 


k 2 ~ h F \ \\ f n +~k x \ G n + ^h > {h„ + 2 W 1 

H n+\ m \ 

1 

H„+- m. 


1 


l 2 =hF 2 \\f n +-k x , G n +-h 


1 


1 


m 1 = hF i /„ +-k x ,\G„ +-li 


1 


Vn 


n„ 


In a similar way k 3 , l 3 , m 3 and k 4 , l 4 , m 4 are calculated following standard 
formulae for the Runge Kutta integration. For example, k 3 is given by 


k 3 - hF i \\ f n + — k 2 U G n + — 1 2 , H n + — m 2 , \t]„ + — 


The functions F u F 2 and F 3 are G, H, ////2 respectively. Then at a distance A T? 
from the wall, we have 


/(At?) =m + G( 0) At? 

(9.26a) 

G(Aij) = G(0) + H( 0) At? 

(9.26b) 

//(At?) = //(0) + /f(0) At? 

(9.26c) 

/T(At?) = - j /(At?) H (At?) 

(9.26d) 


As it has been mentioned earlier /"(0) = H( 0) = A is unknown. It must be 
determined such that the condition = G(°°) = 1 is satisfied. The condition at 
infinity is usually approximated at a finite value of ?7 (around i] = 10). The process 
of obtaining A accurately involves iteration and may be calculated using the 
procedure described below. 

For this purpose, consider Fig. 9.2(a) where the solutions of G versus rj for 
two different values of H( 0) are plotted. The values of G(°°) are estimated from 
the G curves and are plotted in Fig. 9.2(b). The value of H( 0) now can be 
calculated by finding the value H(0 ) at which the line 1-2 crosses the line 
G(oo) = 1. By using similar traingles, it can be said that 


H(0)-H(0) x _ H( 0) 2 -7/(0)! 

1 — G(°°)j G(°°) 2 — G(°°)i 

By solving this, we get H (0). Next we repeat the same calculation as above by 
using H (0) and the better of the two initial values of //(()). Thus we get another 
improvedvalue H (0). This process may continue, that is, we useH (0) and H (0) 
as a pair of values to find more improved values for H( 0), and soforth. It should 
be always kept in mind that for each value of//(0), the curve G( rj) versus rj is to 
be examined to get the proper value of G(°°). 

The functions f{vj), f'(ij ) = G and /"( tj) = H are plotted in Fig. 9.3. The 
velocity components, u and v inside the boundary layer can be computed from 




















The McGraw-Hill Companies 


368 


Introduction to Fluid Mechanics and Fluid Machines 



(a) 



Fig. 9.2 Correcting the initial guess for H(O) 
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Fig. 9.3 f, G, and H distribution in the boundary layer 
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Eqs (9.21a) and (9.21b) respectively. Measurements to test the accuracy of 
theoretical results were earned out by many scientists. In his experiments, J. 
Nikuradse, found excellent agreement with the theoretical results with respect to 
velocity distribution (u/U^) within the boundary layer of a stream of air on a flat 
plate. However, some values of the velocity profile shape f'( Tj ) = u/U,„ = G and 
f"(j ]) = H are given in Table 9.1. 


Table 9.1 

Blasius Velocity Profile G = 

u/U x , f and 

H after Schlichting [2] 

1 

/ 

G 

H 

0 

0 

0 

0.33206 

0.2 

0.00664 

0.06641 

0.33199 

0.4 

0.02656 

0.13277 

0.33147 

0.8 

0.10611 

0.26471 

0.32739 

1.2 

0.23795 

0.39378 

0.31659 

1.6 

0.42032 

0.51676 

0.29667 

2.0 

0.65003 

0.62977 

0.26675 

2.4 

0.92230 

0.72899 

0.22809 

2.8 

1.23099 

0.81152 

0.18401 

3.2 

1.56911 

0.87609 

0.13913 

3.6 

1.92954 

0.92333 

0.09809 

4.0 

2.30576 

0.95552 

0.06424 

4.4 

2.69238 

0.97587 

0.03897 

4.8 

3.08534 

0.98779 

0.02187 

5.0 

3.28329 

0.99155 

0.01591 

8.8 

7.07923 

1.00000 

0.00000 


9.4 WALL SHEAR AND BOUNDARY LAYER 
THICKNESS 


With the profile known, wall shear can be evaluated as 

du 


Tv = A* 


dy 


y = o 


or 


or 


H U °o 



dri 

dy 


r\ = 0 


u x 0.33206 x - 

J(vx)/U M 

[/"(0) = 0.33206 from Table 9.1] 
0.332 p Uj 

7 


or 


(9.27a) 
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and the local skin friction coefficient is 

r* — ^ w 

-pul 

2 

Substituting from (9.27a) we get 
^ _ 0.664 


(9.27b) 


In 1951, Liepmann and Dhawan [3], measured the shearing stress on a flat 
plate directly. Their results showed a striking confirmation of Eq. (9.27). 

Total frictional force per unit width for the plate of length L is 


or 


F= \t w dx 
o 

f 0.332 pul dx 

F ~ r 


or 


or 


F = 


0.332 pU„ x 
x 


yfUjV 


1/2 

1 

2 


Jo 


F= 0.664 x pul 


vL 


l<. 


and the average skin friction coefficient is 

- F 1.328 

C / “ 


Upuh) 


(9.28) 


(9.29) 


where, Rcy = U,„ L tv. 

Since u!U x approaches 1.0 asy —> °<>, it is customary to select the boundary layer 
thickness 8 as that point where u/U rja approaches 0.99. From Table 9.1, u/U m 
reaches 0.99 at rj = 5.0 and we can write 


or 



5.0 x 


(9.30) 


However, the aforesaid definition of boundary layer thickness is somewhat 
arbitrary, a physically more meaningful measure of boundary layer estimation is 
expressed through displacement thickness. 
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Fig. 9.4 (a) Displacement thickness (b) Momentum thickness 

Displacement thickness (8*) is defined as the distance by which the external 
potential flow is displaced outwards due to the decrease in velocity in the bound¬ 
ary layer. 


U^S* = J (LL - u) dy 


Therefore, 


*Whr 


dy 


(9.31) 


Substituting the values of {ulU„) and 77 from Eqs (9.21a) and (9.19) into 
Eq. (9.31), we obtain 


<5 = 


f (1 -/0 dt? - lim [77 - /( 77 )] 

' c/ JJ->“ 


or 


8 = 1.7208 


vx 1.7208 x 


(9.32) 


Following the analogy of the displacement thickness, a momentum thickness 
may be defined. Momentum thickness (8 ) is defined as the loss of momentum 
in the boundary layer as compared with that of potential flow. Thus 


p U* 8*' = j pu ([/„ - u) dy 


or 


j* a 

~ * ul 


1 - 




dy 


(9.33) 


0 “ “ v w “ / 

With the substitution of (u/UJ) and r) from Eq. (9.21a) and (9.19), we can 
evaluate numerically the value of 8 ** for a flat plate as 


or 


o = 


8 = 0.664 


J7'(i-/')dt? 

0.664x 


vx 

U 


VRe 


(9.34) 


The relationships between 8, 8* and <5** have been shown in Fig. 9.4. 
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9.5 MOMENTUM-INTEGRAL EQUATIONS FOR 
BOUNDARY LAYER 


If we are to employ boundary layer concepts in real engineering designs, we need 
to devise approximate methods that would quickly lead to an answer even if the 
accuracy is somewhat less. Karman and Pohlhausen devised a simplified method 
by satisfying only the boundary conditions of the boundary layer flow rather than 
satisfying Prandtl’s differential equations for each and every particle within the 
boundary layer. We shall discuss this method herein. 

Consider the case of steady, two-dimensional and incompressible flow, i.e. we 
shall refer to Eqs (9.10) to (9.14). Upon integrating the dimensional form of 
Eq. (9.10) with respect to y = 0 (wall) to y = 8 (where 8 signifies the interface of 
the free stream and the boundary layer), we obtain 


du du 


dv = J 


1 dp d u ^ 
-- +v 


/ 0 

5 <5 8 

du 


p dx dy 


dy 


or 




du 

V ” 


0 0 - 0 ^ 0 

The second term of the left hand side can be expanded as 


\ v j U d .V = Mo 

0 d y 0 d y 

r du , f <9u | . cm 

\v—ay = v d + \u — ay since — = - 

J Wi? •'Wv l rir 


d\ 


or 


or 


dy 

du 

dy 


du dv 


dx dy 


r du r <9u r cm , 


<9u 

r- 

dx 


Substituting Eq. (9.36) in Eq. (9.35) we obtain 




• du 


1 dp 


du 


dx 


dx 


p dx 


dy 


y =0 


dp 


(9.35) 


(9.36) 


(9.37) 


Substituting the relation between — and the free stream velocity for the 

dr 

inviscid zone in Eq. (9.37) we get 


J OX J Wv J 


du 

dx 


or 




du 

dx 


d £/, 

dx 


d U m 
dr 


dy: 


du 

% 


y =0 


V 


dv = -^ 
P 
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J UZ3 


which is reduced to 


8 i jj 8 

J {[» (U x - u)\ dy} + - - j (U^ - u ) d y = 


P 


o dx LL ” o 

Since the integrals vanish outside the boundary layer, we are allowed to put 

8 =° o. 


J ^ [w (£/ m - w)]dy + J (U «, - u) dy 


or 


d 

dx 


<9x L °° J ' dx 

U 

J [« - i/)]dy + J (C/„ - u) dy = 


P 

P 


Substituting Eq. (9.31) and (9.33) in Eq. (9.38) we obtain 
d \ tt 2 


— \U18' 
dx L 


J dx p 


(9.38) 


(9.39) 


Equation (9.39) is known as momentum integral equation for two dimensional 
incompressible laminar boundary layer. The same remains valid for turbulent 
boundary layers as well. Needless to say, the wall shear stress (z w ) will be 

d£7 


different for laminar and turbulent flows. The term jj 


dx 


signifies 


spacewise acceleration of the free stream. Existence of this term means the 

presence of free stream pressure gradient in the flow direction. For example, we 
d U 


get finite value of U, 


dx 


outside the boundary layer in the entrance region 

dt/. 


of a pipe or a channel. For external flows, the existence of 

on the shape of the body. During the flow over a flat plate, yJ, 
the momentum integral equation is reduced to 

— \ul ,5**1 = — 

dx L J P 


dx 

dLL 

dx 


depends 
= 0 and 

(9.40) 


9.6 SEPARATION OF BOUNDARY LAYER 


It has been observed that the flow is reversed at the vicinity of the wall under 
certain conditions. The phenomenon is termed as separation of boundary layer. 
Separation takes place due to excessive momentum loss near the wall in a 
boundary layer trying to move downstream against increasingpressure, i.e., dp/dx 
> 0, which is called adverse pressure gradient. Figure 9.5 shows the flow past a 
circular cylinder, in an infinite medium. Up to 9 = 90°, the flow area is like a 
constricted passage and the flow behaviour is like that of a nozzle. Beyond 
9 = 90° the flow area is diverged, therefore, the flow behaviour is much similar to 
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a diffuser. This dictates the inviscid pressure distribution on the cylinder which is 
shown by a firm line in Fig. 9.5. Flere p„, and U m are the pressure and velocity in 
the free stream and p is the local pressure on the cylinder. 

Boundary 



6 ->- 

0° 45° 90° 135° 180° 



Fig. 9.5 Flow separation and formation of wake behind a circular cylinder 

Consider the forces in the flow field. It is evident that in the inviscid region, the 
pressure force and the force due to streamwise acceleration are acting in the same 
direction (pressure gradient being negative/favourable) until 9 = 90°. Beyond 
9= 90°, the pressure gradient is positive or adverse. Due to the adverse pressure 
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gradient the pressure force and the force due to acceleration will be opposing each 
other in the inviscid zone of this part. So long as no viscous effect is considered, 
the situation does not cause any sensation. However, in the viscid region (near the 
solid boundary), up to 9 = 90°, the viscous force opposes the combined pressure 
force and the force due to acceleration. Fluid particles overcome this viscous 
resistance. Beyond 9 = 90°, within the viscous zone, the flow structure becomes 
different. It is seen that the force due to acceleration is opposed by both the 
viscous force and pressure force. Depending upon the magnitude of adverse 
pressure gradient, somewhere around 6 = 90°, the fluid particles, in the boundary 
layer are separated from the wall and driven in the upstream direction. However, 
the far field external stream pushes back these separated layers together with it 
and develops a broad pulsating wake behind the cylinder. Now let us look at the 
mathematical explanation of flow-separation. Following the foregoing 
observation, the point of separation may be defined as the limit between forward 
and reverse flow in the layer very close to the wall, i.e., at the point of separation 

f du^ 


dy 


= 0 


(9.41) 


This means that the shear stress at the wall, T w = 0. But at this point, the adverse 
pressure continues to exist and at the downstream of this point the flow acts in a 
reverse direction resulting in a back flow. We can also explain flow separation 
using the argument about the second derivative of velocity u at the wall. From the 
dimensional form of the momentum Eq. (9.10) at the wall, where u = v = 0, we 
can write 


d v 


dy 2 


/y = 0 


dp 

dx 


(9.42) 


dp 


Consider the situation due to a favourable pressure gradient were — < 0. 

dx 

From Eq. (9.42) we have, (d 2 u/dy 2 ) wall < 0. As we proceed towards the free 


stream, the velocity u approaches U,„ aysmptotically, so dn/dy decreases at a 
continuously lesser rate in y direction. This means that (c^u/dy 2 ) remains less 
than zero near the edge of the boundary layer. Finally it can be said that for a 
decreasing pressure gradient, the curvature of a velocity profile (c^u/dy 2 ) is 
always negative as shown in (Fig. 9.6a). Next consider the case of adverse 
pressure gradient, dp/dx > 0. From Eq. (9.42), we observe that at the boundary, 
the curvature of the profile must be positive (since dp/chc > 0). 

However, near the interface of boundary layer and free stream the previous 
argument regarding duldy and d 2 u/dy 2 still holds good and the curvature is 
negative. Thus we observe that for an adverse pressure gradient, there must exist 
a point for which d 2 u/dy~ = 0. This point is known as point of inflection of the 
velocity profile in the boundary layer as shown in Fig. 9.6b. However, point of 
separation means dn/dy = 0 at the wall. In addition, Eq. (9.42) depict d 2 uld\’ 2 > 0 
at the wall since separation can only occur due to adverse pressure gradient. But 
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we have already seen that at the edge of the boundary layer, d 2 u/dy 2 < 0. It is 
therefore, clear that if there is a point of separation, there must exist a point of 
inflection in the velocity profile. 



(a) 




Fig. 9.6 Velocity distribution within a boundary layer 

(a) Favourable pressure gradient, — < 0 

dx 

(b) adverse pressure gradient, — > 0 

dx 

Let us reconsider the flow past a circular cylinder and continue our discussion 
on the wake behind a cylinder. The pressure distribution which was shown by the 
firm line in Fig. 9.5 is obtained from the potential flow theory. However, some¬ 
where near 0 = 90° (in experiments it has been observed to be at 0 = 81°), the 
boundary layer detaches itself from the wall. Meanwhile, pressure in the wake 
remains close to separation-point-pressure since the eddies (formed as a conse¬ 
quence of the retarded layers being carried together with the upper layer through 
the action of shear) cannot convert rotational kinetic energy into pressure head. 
The actual pressure distribution is shown by the dotted line in Fig. 9.5. Since 
the wake zone pressure is less than that of the forward stagnation point (pres¬ 
sure at point A in Fig. 9.5), the cylinder experiences a drag force which is 
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basically attributed to the pressure difference. The drag force, brought about by 
the pressure difference is known as form drag whereas the shear stress at the wall 
gives rise to skin friction drag. Generally, these two drag forces together are 
responsible for resultant drag on a body. 

9.7 KARMAN-POHLHAUSEN APPROXIMATE 
METHOD FOR SOLUTION OF 
MOMENTUM INTEGRAL EQUATION 
OVER A FLAT PLATE 


The basic equation for this method is obtained by integrating the x direction 
momentum equation (boundary layer momentum equation) with respect toy from 
the wall (at y = 0 ) to a distance <5 (x) which is assumed to be outside the boundary 
layer. With this notation, we can rewrite the Karman momentum integral equation 
(9.39) as 

tt 2 d<5 /- s .** s*\tt K 

u -~i r + ( M (943) 

The effect of pressure gradient is described by the second term on the left hand 
side. For pressure gradient surfaces in external flow or for the developing sections 
in internal flow, this term will be retained and will contribute to the pressure 
gradient. However, we assume a velocity profile which is a polynomial of 1 ) = 
y/8. As it has been seen earlier, t] is a form of similarity variable. This implies 
that with the growth of boundary layer as distance x varies from the leading edge, 
the velocity profile (u/U^) remains geometrically similar. We choose a velocity 
profile in the form 


= a 0 + a, t) + a 2 r\ 2 + a 3 rf (9.44) 

In order to determine the constants a 0 , a l9 a 2 and a 3 we shall prescribe the 
following boundary conditions 


at 


at 


at 


at 


y = o 


y = 0 , u = 0 
d 2 ~ 


= 0 


dy 2 

y = 8, u = U x 


or at ri = 0 , -= 0 


or at ri = 0, - 9 (u/U^ ) = 0 

drj- 


or 


x du n 

y = o, — = 0 

dy 


or 


at r\ = 1, 
at 77 = 1 , 


= 1 


= 0 


u, 

d{u!U „) 
dt] 

These requirements will yield 

a 0 = 0, a 2 = 0, aj + 3a 3 = 0 and a x + a 3 = 1 
Finally, we obtain the following values for the coefficients in Eq. (9.44), 


(9.45a) 

(9.45b) 

(9.45c) 

(9.45d) 
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3 1 

a 0 = 0, a x = a 2 = 0 and a 3 = - — 

and the velocity profile becomes 

U 3 1 3 

-= —rj -rf 

U 2 2 


For flow over a flat plate, —— = 0, hence U,„ —^ 

dr dr 

Eq. (9.43) reduces to 

dS ** r w 


dx put 

Again from Eq. (9.33), the momentum thickness is 


(9.46) 

= 0 and the governing 

(9.47) 


or 


or 


or 


r u 

o = - 

J 0 U ~ 
8 

s"-\ 


u 

uL 


1 -— 

t/ 


u,. 


d.V 


d.V 




3 1 3^ 

r? - rj 

2 2 


dr] 


§** = s 
280 


The wall shear stress is given by 

du 




dy 


y = o 


or 


or 


Tv = M 


T„, = 


8 drj 

28 


u A r] -\ r? 


Jrj=0 


Substituting the values of <5** and z w in Eq. (9.47) we get, 
39 d 8 2>p U „ 


or 


280 dr 28pUl 
■140 p 




or 


13 pU a 
140 vx 


dr F (. 


C, 


8^_ _ 

2 13 U„ 

where Q is any arbitrary unknown constant. 


(9.48) 
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The condition at the leading edge (at x = 0, <5 = 0) yields 

Cj = 0 

Finally we obtain, 

^ = 280 vx 
13 C/ 


(9.49) 


or 


or 


<5=4.64 


vx 


<5 = 


4.64x 

VRe7 


(9.50) 


This is the value of boundary layer thickness on a flat plate. Although, the 
method is an approximate one, the result is found to be reasonably accurate. The 
value is slightly lower than the exact solution of laminar flow over a flat plate 
given by Eq. (9.30). As such, the accuracy depends on the order of the velocity 
profile. It maybe mentioned that instead of Eq. (9.44), we can as well use a fourth 
order polynomial as 


-jj— = a 0 + a x ri + a 2 rj 2 + a 3 rj 3 + a 4 ;j 4 


(9.51) 


In addition to the boundary conditions in Eq. (9.45), we shall require another 
boundary condition 


at 


y= 8, 


d 2 


u , d'-{u/U oa ) 

Y — 0 or at 77 =1 , - - = 0 


dy l ‘ drj 2 

To determine the constants as a 0 = 0, a l = 2, a 2 = 0, a 3 = -2 and a 4 = 1. Finally the 
velocity profile will be 

— = 2 77 — 277 s + t/ 4 

Subsequently, for a fourth order profile the growth of boundary layer is given by 
<5=4^ (9-52) 

This is also very close to the value of the exact solution. 

9.8 INTEGRAL METHOD FOR NON-ZERO 
PRESSURE GRADIENT FLOWS 


A wide variety of “integral methods” in this category have been discussed by 
Rosenhead [4]. The Thwaites method [5] is found to be a very elegant method, 
which is an extension of the method due to Holstein and Bohlen [ 6 ]. We shall 
discuss the Holstein-Bohlen method in this section. 

This is an approximate method for solving boundary layer equations for two- 
dimensional generalized flow. The integrated Eq. (9.39) for laminar flow with 
pressure gradient can be written as 
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/[t/ 2 <T]+<rc4^ & — 

dx L J dx p 


or 


U‘ 


d5 

dx 


(2 S**+S*)U^ = ^- 

J dx p 


(9.53) 


The velocity profile the boundary layer is considered to be a fourth-order 
polynomial in terms of the dimensionless distance i) = y!5, and is expressed as 

u/U = ap + bp 2 + cp 2 + dp 4 
The boundary conditions are 


p = 0 : u = 0 , v = 0 


p = 1: u = U, 


v d 2 u _ id p _ dU 
e2 “ ~n 7TV “ U dx 


5 2 dp 2 P dx 


dii _ „ d 2 u 
dp ’ dp 2 


= 0 


A dimensionless quantity, known as shape factor is introduced as 

x= s ±w 

V dx 

The following relations are obtained 

A A A 

a = 2 + ~, b= -, c = -2 H—, d= 1 - 

6 2 2 

Now, the velocity profile can be expressed as 

u/U = F{p) + XG(p), 

where 

F( p) = 2 p - 2 p 2 + p 4 , G( p) = i r?( 1 - p) 2 

6 

The shear stress T w = fiidu/dy)^ is given by 
pU 6 

We use the following dimensionless parameters, 

o** 

T,.,S 


L = 


pU 


(8*y du 

K ~~V~ d7 

The integrated momentum Eq. (9.53) reduces to 

d<T 


U- 


dx 


(2 +H)^ = ^ 

dx o 


(9.54) 


A 


(9.55) 


(9.56) 

(9.57) 

(9.58) 

(9.59) 
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or 




2[L - K (H + 2)] 


(9.60) 


The parameter L is related to the skin friction and K is linked to the pressure 
gradient. If we take K as the independent variable, L and H can be shown to be the 
functions of K since 


<5 


i 

J [l-F(t 7 )-AG(t 7 )]dt 7 


10 


A 

120 


= —- — ( 9 - 61 ) 


l 

J (Fin) + m-n) (i - F(n) - M(n))dn 


K = 


o 

37 A 
315 945 

[<5**] 2 


A 


9072 


;2 


-A = A 


37 

315 


A 

945 


5 


9072 


(9.62) 


(9.63) 


Therefore, 


L = \2 + — 


6)8 


-4 


37 

315 


A 

945 


A* 


2 A 


// - 


(3/10) — (A/120) 


9072^ 
=/ 2 (^) 


=/,m 


<5" (37/315) - (A/945) - (A 2 /9072) 

The right-hand side of Eq. (9.60) is thus a function of K alone. Walz [7] pointed 
out that this function can be approximated with a good degree of accuracy by a 
linear function of K so that 

2 [L - K(H +2 )\=a-bK 
Equation (9.60) can now be written as 

/^rs-i2A [/[S’*] 2 1 d U 


_d 

dx 


U[8 f 


= a 


(b- 1 ) 


U dx 


V v 

Solution of this differential equation for the dependent variable (£/[<5**] 2 /v) 
subject to the boundary condition (7=0 when x = 0, gives 


£/[<n 2 


X 



0 

With a = 0.47 and 6 = 6 , the approximation is particularly close between the 
stagnation point and the point of maximum velocity. Finally the value of the 
dependent variable is 
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= V s d* (9-64) 

0 

By taking the limit of Eq. (9.64), according to L’Hospital’s rule, it can be shown 
that 

[<5**] 2 | x = o = 0-47 v/6U'(0) 

This corresponds to K = 0.0783. It may be mentioned that [<5 ] is not equal to 
zero at the stagnation point. If ([<5**] 2 /v) is determined from Eq. (9.64), K{x) can 
be obtained from Eq. (9.58). Table 9.2 gives the necessary parameters for 
obtaining results, such as velocity profile and shear stress T w . The approximate 
method can be applied successfully to a wide range of problems. 


Table 9.2 Auxiliary functions after Holstein and Bohlen [6] 


A 

K 

fM 


12 

0.0948 

2.250 

0.356 

10 

0.0919 

2.260 

0.351 

8 

0.0831 

2.289 

0.340 

7.6 

0.0807 

2.297 

0.337 

7.2 

0.0781 

2.305 

0.333 

7.0 

0.0767 

2.309 

0.331 

6.6 

0.0737 

2.318 

0.328 

6.2 

0.0706 

2.328 

0.324 

5.0 

0.0599 

2.361 

0.310 

3.0 

0.0385 

2.427 

0.283 

1.0 

0.0135 

2.508 

0.252 

0 

0 

2.554 

0.235 

- 1 

-0.0140 

2.604 

0.217 

-3 

- 0.0429 

2.716 

0.179 

-5 

-0.0720 

2.847 

0.140 

-7 

-0.0999 

2.999 

0.100 

-9 

-0.1254 

3.176 

0.059 

- 11 

-0.1474 

3.383 

0.019 

- 12 

-0.1567 

3.500 

0 


As mentioned earlier, K and A are related to the pressure gradient and the shape 
factor. Introduction of K and A in the integral analysis enables extension of 
Kannan-Pohlhausen method for solving flows over curved geometry. However, 
the analysis is not valid for the geometries, where A < - 12 and A > + 12. 

9.9 ENTRY FLOW IN A DUCT 


Growth of boundary layer has a remarkable influence on flow through a constant 
area duct or pipe. Let us consider a flow entering a pipe with uniform velocity. 
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The boundary layer starts growing on the wall at the entrance of the pipe. 
Gradually it becomes thicker in the downstream and the flow becomes fully 
developed when the boundary layers from the wall meet at the axis of the pipe 
(Fig. 9.7). The velocity profile is nearly rectangular at the entrance and it 
gradually changes to a parabolic profile at the fully developed region. Before the 
boundary layers from the periphery meet at the axis, there prevails a core region 
which is uninfluenced by viscosity. Since the volume-flow must be same for every 
section and the boundary-layer thickness increases in the flow direction, the 
inviscid core accelerates, and there is a corresponding fall in pressure. It can be 
shown that for laminar incompressible flows, the velocity profile approaches the 
parabolic profile through a distance Le from the entry of the pipe, which is given 
by 

— = 0.05 Re, where Re = 

D v 

For a Reynolds number of 2000, this distance, which is often referred to as the 
entrance length is about 100 pipe-diameters. For turbulent flows, the entrance 
region is shorter, since the turbulent boundary layer grows faster. 

At the entrance region, the velocity gradient is steeper at the wall, causing a 
higher value of shear stress as compared to a developed flow. In addition, 
momentum flux across any section in the entrance region is higher than that 
typically at the inlet due to the change in shape of the velocity profile. Arising out 
of these, an additional pressure drop is brought about at the entrance region as 
compared to the pressure drop in the fully developed region. 


Le 



Fig. 9.7 Development of boundary layer in the entrance region of a duct 

9.10 CONTROL OF BOUNDARY LAYER 
SEPARATION 

It has already been seen that the total drag on a body is attributed to form drag and 
skin friction drag. In some flow configurations, the contribution of form drag 
becomes significant. In order to reduce the form drag, the boundary layer 
separation should be prevented or delayed so that somewhat better pressure 
recovery takes place and the form drag is reduced considerably. There are some 
popular methods for this purpose which are stated as follows. 
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(i) By giving the profile of the body a streamlined shape as shown in 
Fig. 9.8. This has an elongated shape in the rear part to reduce the 
magnitude of the pressure gradient. The optimum contour for a 
streamlined body is the one for which the wake zone is very narrow 
and the form drag is minimum. 



Fig. 9.8 Reduction of drag coefficient (C D ) by giving the profile a streamlined shape 


(ii) The injection of fluid through porous wall can also control the 
boundary layer separation. This is generally accomplished by blowing 
high energy fluid particles tangentially from the location where 
separation would have taken place otherwise. This is shown in 
Fig. 9.9. The injection of fluid promotes turbulence and thereby 
increases skin friction. But the form drag is reduced considerably due 
to suppression of flow separation and this reduction can be of 
significant magnitude so as to ignore the enhanced skin friction drag. 



9.11 MECHANICS OF BOUNDARY LAYER 
TRANSITION 

One of the interesting problems in fluid mechanics is the physical mechanism of 
transition from laminar to turbulent flow. The problem evolves about the 
generation of both steady and unsteady vorticity near a body, its subsequent 
molecular diffusion, its kinematic and dynamic convection and redistribution 
downstream, and the resulting feedback on the velocity and pressure fields near 
the body. We can perhaps realise the complexity of the transition problem by 
examining the behaviour of a real flow past a cylinder. 

Figure 9.10 (a) shows the flow past a cylinder for a very low Reynolds number 
(~1). The flow smoothly divides and reunites around the cylinder. At a Reynolds 
number of about 4, the flow separates in the downstream and the wake is formed 
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by two symmetric eddies. The eddies remain steady and symmetrical but grow in 
size up to a Reynolds number of about 40 as shown in Fig. 9.10(b). 

When the Reynolds number crosses 40, oscillation in the wake induces 
asymmetry and finally the wake starts shedding vortices into the stream. This 
situation is termed as onset of periodicity as shown in Fig. 9.10(c) and the wake 
keeps on undulating up to a Reynolds number of 90. As the Reynolds number 
further increases, the eddies are shed alternately from a top and bottom of the 
cylinder and the regular pattern of alternately shed clockwise and counter¬ 
clockwise vortices form Von Karman vortex street as in Fig. 9.10(d). However, 
periodicity is eventually induced in the flow field with the vortex-shedding 
phenomenon. The periodicity is characterised by the frequency of vortex 
shedding /. In non-dimensional form, the vortex shedding frequency is expressed 
as fD/U ve p known as the Strouhal number named after V. Strouhal, a German 
physicist who experimented with wires singing in the wind. The Strouhal number 
shows a slight but continuous variation with Reynolds number around a value of 
0.21. At about Re = 500, multiple frequencies start showing up and the wake 
tends to become turbulent. 


Re ~ 1 


Re = U^D/V 

Re ~ 4 to about 40 



Fig. 9.10 Influence of Reynolds number on wake-zone aerodynamics 

An understanding of the transitional flow processes will help in practical 
problems either by improving procedures for predicting positions or for 
determining methods of advancing or retarding the transition position. 

The critical value at which the transition occurs in pipe flow is Re cr = 2300. 
The actual value depends upon the disturbance in flow. Some experiments have 
shown the critical Reynolds number to reach as high as 40,000. The precise upper 
bound is not known, but the lower bound appears to be Re„. = 2300. Below this 
value, the flow remains laminar even when subjected to strong disturbances. 
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For 2300 < Rc rr < 2600, the flow alternates randomly between laminar and 
partially turbulent. Near the centerline, the flow is more laminar than turbulent, 
whereas near the wall, the flow is more turbulent than laminar. For flow over a 
flat plate, turbulent regime is observed between Reynolds numbers ([/«, x/v) of 
3.5 x 10 5 and 10 6 . 

9.12 SEVERAL EVENTS OF TRANSITION 

Transitional flow consists of several events as shown in Fig. 9.11. Let us consider 
the events one after another. 

1. Region of instability of small wavy disturbances Consider a laminar flow over 
a flat plate aligned with the flow direction (Fig. 9.11). It has been seen that in the 
presence of an adverse pressure gradient, at a high Reynolds number (water 
velocity approximately 9-cm/sec), two-dimensional waves appear. These waves 
can be made visible by a method known as tellurium method. In 1929, Tollmien 
and Schlichting predicted that the waves would form and grow in the boundary 
layer. These waves are called Tollmien-Schlichting wave. 

2. Three-dimensional waves and vortex formation Disturbances in the free 
stream or oscillations in the upstream boundary layer can generate wave growth, 
which has a variation in the spanwise direction. This leads an initially two- 
dimensional wave to a three-dimensional form. In many such transitional flows, 
periodicity is observed in the spanwise direction. This is accompanied by the 
appearance of vortices whose axes lie in the direction of flow. 

5. Peak-Valley development with streamwise vortices As the three-dimensional 
wave propagates downstream, the boundary layer flow develops into a complex 
streamwise vortex system. However, within this vortex system, at some spanwise 
location, the velocities fluctuate violently. These locations are called peaks and 
the neighbouring locations of the peaks are valleys (Fig. 9.12). 

4. Vorticity concentration and shear layer development At the spanwise 
locations corresponding to the peak, the instantaneous streamwise velocity 
profiles demonstrate the following. Often, an inflexion is observed on the velocity 
profile. The inflectional profile appears and disappears once after each cycle of 
the basic wave. 

5. Breakdown The instantaneous velocity profiles produce high shear in the 
outer region of the boundary layer. The velocity fluctuations develop from the 
shear layer at a higher frequency than that of the basic wave. However, these 
velocity fluctuations have a strong ability to amplify any slight three- 
dimensionality, which is already present in the flow field. As a result, a staggered 
vortex pattern evolves with the streamwise wavelength twice the wavelength of 
Tollmien-Schlichting wavelength. The spanwise wavelength of these structures 
is about one-half of the streamwise value. This is known as breakdown. Klebanoff 
et al. [8] refer to the high frequency fluctuations as hairpin eddies. 
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6. Turbulent-spot development The hairpin-eddies travel at a speed grater than 
that of the basic (primary) waves. As they travel downstream, eddies spread in 
the spanwise direction and towards the wall. The vortices begin a cascading 
breakdown into smaller vortices. In such a fluctuating state, intense local changes 
occur at random locations in the shear layer near the wall in the form of turbulent 
spots. Each spot grows almost linearly with the downstream distance. The 
creation of spots is considered as the main event of transition. 



Fig. 9.12 Cross-stream view of the streamwise vortex system 


Summary 

• The boundary layer is the thin layer of fluid adjacent to the solid 
surface. Phenomenologically, the effect of viscosity is very prominent 
within this layer. 

• The main-stream velocity undergoes a change from zero at the solid- 
surface to the full magnitude through the boundary layer. 

• Effectively, the boundary layer theory is a complement to the inviscid 
flow theory. 

• The governing equation for the boundary layer canbe obtained through 
correct reduction of the Navier-Stokes equations within the thin layer 
referred above. There is no variation in pressure in y direction within 
the boundary layer. 

• The pressure is impressed on the boundary layer by the outer inviscid 
flow which can be calculated using Bernoulli’s equation. 

• The boundary layer equation is a second order non-linear partial 
differential equation. The exact solution of this equation is known as 
similarity solution. For the flow over a flat plate, the similarity 
solution is often referred to as Blasius solution. Complete analytical 
treatment of this solution is beyond the scope of this text. However, the 
momentum integral equation can be derived from the boundary layer 
equation which is amenable to analytical treatment. 
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• The solutions of the momentum integral equation are called approx¬ 
imate solutions of the boundary layer equation. 

• The boundary layer equations are valid up to the point of separation. 
At the point of separation, the flow gets detached from the solid surface 
due to excessive adverse pressure gradient. 

• Beyond the point of separation, the flow reversal produces eddies. 
During flow past bluff-bodies, the desired pressure recovery does not 
take place in a separated flow and the situation gives rise to pressure 
drag or form drag. 
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Solved Examples 


Ex amp le 9.1 _Water flows over a flat plate at a free stream velocity of 0.15 m/s. 

There is no pressure gradient and laminar boundary layer is 6 mm thick. Assume a 
sinusoidal velocity profile 

u _ . k f v 
U„ Sm 2{S 

For the flow conditions stated above, calculate the local wall shear stress and skin friction 
coefficient. 


[p = 1.02 x 10 3 kg/ms, p = 1000 kg/m J ] 
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u» 


Solution 


or 


du uUao d(u/U x ) 

T= u — — — ---- 

dv 8 d{y/8) 

k (n y'\_ 1-57p£/ c 


8 2 


— cos — ~ = 


2 8 


n y 

-cos! 


= T 


1.57 /iU„ 


\y=0 ' 


1.57 x 1.02 x 10“ 3 x 0.15 A n . XT/ 2 
T„, =--—— ^3 -= 0.04 N/m 


cv= 


6x10“ 

2 x 0.04 


/ 1 


pU 2 1000 x (0.15) z 


= 3.5 x 10“ 


Ex amp le 9.2 Air at standard conditions flows over a flat plate. The freestream 
speed is 3 m/sec. Find (5and t w atx = 1 m from the leading edge (assume a cubic velocity 
profile). For air, v= 1.5 X 10~ 5 m 2 /sec and p = 1.23 kg/m 3 . 


Solution Applying the results developed in section 9.7 for cubic velocity profile and 
the growth of boundary layer, we can write 


and 


u 3 i 3 , y t 

- = — 77 - 77 , where 77 = — at any x 

U x 2 2 8 

8 4.64 




For air with v = 1.5 x 10 ‘ m 2 /s, the local Reynolds number at x is 


R e A = 
8 = 


U 0 „x 

v 


4.64 x 

VReJ 


3x1 

1.5 x 10~ 5 
4.64x1 
3/2 x 10 5 


= 2x 10 5 

^m = 0.0103 m= 10.3 mm 


T w = M 


du 

dy 


HU" „ 


y =0 


"3 1 3 

-77-77 

d771_2 2 


T]-0 


2/j.U^ _3pvU m 
28 ~ 28 


= 3 x 1,23 x 1,5 x 10~ 5 x 3 
vv 2x0.0103 

T„, = 8.06 x 10 -3 N/m 2 


Example 9.3_ Air moves over a flat plate with a uniform free stream velocity of 10 
m/s. At a position 15 cm away from the front edge of the plate, what is the boundary layer 
thickness? Use a parabolic profile in the boundary layer. For air, v= 1.5 x 10 5 m 2 /s and 
p= 1.23 kg/m 3 . 

Solution For a parabolic profile let us take 
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= a + by + ct 2 
CL 


The boundary conditions are 


at 

y= 

0 , 

u = 0 

at 

y= 

8, 

u = CL 



8, 

du 

at 

y= 

~dv = 


Evaluating the constants we get 


CL 


Now 


*w= P 


8 

du 

dy 


7L=2 t “ T =2f|-»r 


/t(/„ d(u!UJ) 


y =o 




7] — 0 


pCL d(277 — 77 ) 


5 d/7 

Applying momentum integral equation 

1 


2pCL 


7] — 0 


rr2 do f w L w 


2pC/ c 


dx J CL 
0 

d<5 ' 


dp 


s pUi — J ( 2/7 -/ 7 2 )(1 - 2/7 + ? 7 2 ) d /7 


or 


2pCL = d5 

ctpC/i dx 


1 

J (277 - 5 / 7 2 + 4 / 7 3 - p 4 ) dr) 


2 p _ 2 dc> 
15 d7 


<5pL 


<5d5 = dx 


pU„ 

8 2 _ 15p 


x + C 


2 pCL 

It is assumed that at x = 0, 8=0 which yields C = 0. Thus 


8 = 


or 


30 p 

ptf- ' 

30 p 


5.48 


pC/ M x 


In this problem, 


_ 10 x 15 x 10 ~ 2 . ,. 5 

Re =- f— = 1 x 10 

1.5 x 10 ~~ 5 




























The McGraw-Hill Companies 




Laminar Boundary Layers 


391 


or 


5.48 

8= , - x 15 cm = 0.259 cm 

8= 2.59 mm 


Ex amp le 9 .4 Air moves over a 10 m long flat plate. The transition from laminar to 
turbulent flow takes place between Reynolds numbers of 2.5 x 10 6 and 3.6 x 10 6 . What 
are the minimum and maximum distance from the front edge of the plate along which one 
expect laminar flow in the boundary layer? The free stream velocity is 30 m/s and 


v= 1.5 x 10 5 m 2 /s. 


Solution We can see that the range of Reynolds numbers for laminar flow is 2.5 x 10 6 
to 3.6 x 10 6 


For the lower limit, 


2.15x10 = 


30x 


or x n 

for the upper limit, 


1.5 x 10~ 
, = 1.075 m 


3.6 x 10~ 6 = 


30x 


or 


1.5 x 10~ 
= 1.8m 


Exam ple 9 .5 Water at 15 °C flows over a flat plate at a speed of 1.2 m/s. The plate 
is 0.3 m long and 2 m wide. The boundary layer on each surface of the plane is laminar. 

Assume the velocity profile is approximated by a linear expression for which — = ? . 

x jRe x 

Determine the drag force on the plate. For water v= 1.1 x 10" 6 m 2 /s, p = 1000 kg/m 3 . 
Solution On a flat plate, the drag is due to skin friction acting on each side of the plate 

L 


Fd = 2 


T ux b dx 


For linear profile = — and T„, = u — 

U x 8 w ^ dy 


y =o 


r = pU„ d(U/U x ) 


8 

L 


d(y!8) 


F d = 2 




y=0 

L 


b dx = 2 


8 


o 

2 juU„ 
3.46 


3.46 V vx 


■b dx 


L 


dx 
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Therefore, 

Thus, 


2 nUJ> 
3.46 



4 nUJb LL_L _ ytUJ) 
3.46 V v 3.46 


VRe7 


Re, = 


UL 


L2X03 = 3.27 x 10 5 


(Re z ) 1/2 = 572 


1.1x10 


4 x 1.1 x 10~ 3 x 1.2x2x572 
3.46 


-6 


F d = 


= 1.745 N 


Ex amp le 9.6 _ Air is flowing over a thin flat plate which is 1 m long and 0.3 m wide. 
At the leading edge, the flow is assumed to be uniform and U„ = 30 m/s. The flow 
condition is independent of z (see Fig. 9.13). Using the control volume abed, calculate 
the mass flow rate across the surface ab. Determine the magnitude and direction of the x 
component of force required to hold the plate stationary. The velocity profile at be is 
given by 



and 8=4 mm. Density of air = 1.23 kg/m 3 and V= 1.5 x 10 5 m 2 /sec. 



► 

x 


Fig. 9.13 Control volume of interest on flat plate 



U 1 

Solution At be, jj- = 2 q q ; i] = y/8 
Steady state continuity equation is given by 

J pV-dA =0 

s 


S 

-pU ^ bd + J pit b dy + m ab = 0 
o 


or 
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but 


S 1 

J pu b dy = pU aa b8 J (2 77 - rj 2 )dri = pU^ b8 


2 rf 

71 "T 


-it 


= -pU x b8 
3 


m ab = 


= pU^bd-- pU^ bS = ^pU oa b8 


or m ab = — x 1.23x30x0.3x0.004 = 0.0147 kg/s 

Steady state momentum equation is given by 

J upV-dA = F sx 


cs 


or 

£ 

js 

II 

X 

, b8 } + u t 

But, 

u da = u ab = U„, 

and 


S 1 

J u pu b dy = pU 2 „ b8 J ( 2 rj - rffdT) 

0 

= puibs 


4 2 4 77 

- n- f] + — 
3 3 


= - pUlbS 
15 F 


Thus, 


15 


Hence, 


R x = -pU 2 „ b8+ y put b8+ ^ pUi b8 = - ^ pU'i b8 

R x = -y {1.23 x(3) 2 x 0.3x0.004} 

= - 0.177 N (the force must be applied to the CV by 
the plate) 

F = R = 0 1 77 N 7tn the leftl 
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Exercises 


9.1 Choose the correct answer 

(i) The laminar boundary layer thickness on a flat plate varies as 

(a) (b) x^ (c) x (1/2 > (d) x 2 

(ii) The turbulent boundary layer thickness on a flat plate varies as 

(a) x< +1/2 > (b) x (4/5) (c) x< 1/7 > (d) 

(iii) The injection of air through a porous wall can control the boundary layer 
separation on the upper curved surface of an aerofoil wing. The injection of 
fluid also promotes turbulence. The final result is 

(a) increase in the skin friction and decrease in the form drag 

(b) increase in the form drag and decrease in the skin friction 

(c) decrease in both the skin friction and form drag 

(d) increase in both the skin friction and form drag 
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(iv) In the entrance region of a pipe, the boundary layer grows and the inviscid 
core accelerates. This is accompanied by a 

(a) rise in pressure 

(b) constant pressure gradient 

(c) fall in pressure in the flow direction 

(d) pressure pulse 

(v) Flow separation is caused by 

(a) reduction of pressure to vapour pressure 

(b) a negative pressure gradient 

(c) a positive pressure gradient 

(d) the boundary layer thickness reducing to zero. 

(vi) At the point of separation 

(a) shear stress is zero 

(b) velocity is negative 

(c) pressure gradient is negative 

(d) shear stress is maximum 

(vii) Choose the expression for the momentum thickness of an incompressible 
boundary layer 


(a) 


5.0 x 

VRe7 



1 -- 


t/„ 


dy 



\ 

U 



dy 


(d) \{u!Uj)dy 
0 


(viii) 


For cross flow over a circular cylinder at a Reynolds number 



U.„D 

v 


greater than 200 , the wake is 

(a) at a pressure equal to the forward stagnation point 

(b) at a pressure lower than the forward stagnation point 

(c) the principal cause of skin friction 

(d) at a pressure higher than the forward stagnation point 

9.2 Nikuradse, a well known student of Prandtl, obtained experimental data for 
laminar flow over a flat plate placed at zero angle of attack. His measurements 
suggest 



where, y is the perpendicular distance from the plate. The local velocity is u. 
Evaluate a and b from physical boundary conditions. Obtain the expressions for 
the boundary layer thickness 8, displacement thickness 8 , momentum thickness 
8 and the shear stress T w on the surface of the plate. 

A ns. (a = 3/2, b = -1/2, 8/x = 4.64 /(Re/' 5 , 
(< 57 * = 1.73 (Re/' 5 , t w = 0.323 (Re/ 5 Ulx 
9.3 Given the choice between cos Ay and sin Ay as velocity profiles, which one would 
you prefer? To determine choice of the profile, find the displacement thickness, 
momentum thickness, wall shear stress and boundary layer thickness, from the 
momentum integral equation for flow over a flat plate. 

Ans. (sin Ay, <5* = 0.363 8, <5" = 0.137 8, x w = iqi UJ28, 8/x = 4.8/(Re/' 5 
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9.4 For the laminar flow over a flat plate, the experiments confirm the velocity profile 

\3 

. For the turbulent flow over a flat plate, the 


u_ _ 3 [ Z|_If Z 

U 2{S 2U 


experimental observations over a range of Reynolds number suggest 


u _ j y 
U ' 


1/7 


Find the ratio of (o IS) for laminar and turbulent cases. 

Ans. (<57<S) Iam = 3/8, (S/8) tm = 1/8 


9.5 Assuming the velocity profile —— = tanh —-— for the boundary layer over a 

U„ a(x) 

flat plate at zero incidence, find the relations for 8, 8* S*" and T w . Check whether 
the profile satisfies all the boundary conditions. 

Note: a(x) ^ 5 (x) where 8 is such that at y = 8,u = 0.99 t/„ 

Ans. (8!x = 6.76/(ReJ 0 ' 5 , Six = 1.77/(Re /' 5 <Tlx = 0.783/(ReJ°' 5 , 

T w = /tt/„(Re x )°' 5 /2.553x 

9.6 An approximate expression for the velocity profile in a steady, 2-D, 
incompressible boundary layer is 

-jj— = l-e~ v +k[\-e~ ri -sin-^J, for 0 < r\ < 3 

= 1 - e n - ke^ n , for 77 > 3 

where 77 = y/8{x). Show that the profile satisfies the following boundary 
conditions 

aty = 0 , m = 0 


9.7 


9.8 


rr dll 
at y = 00, j; = , - 

dy 



= 0 


Also find out k from an appropriate boundary condition. 

Ans. (k= ( S 2 /v) (dUJdx) -1) 
Water of kinematic viscosity v = 1.02 x 10~ 6 m 2 /s is flowing steadily over a 
smooth flat plate at zero angle of attack with a velocity 1.6 m/s. The length of the 
plate is 0.3 m. Calculate (a) the thickness of boundary layer at 15 cm from the 
leading edge (b) the rate-of-growth of boundary layer at 15 cm from the leading 
edge, and (c) the total drag coefficient on one side of the plate. Assume a 
parabolic velocity profile. 


Ans. (8= 1.7 mm, d<S/dx = 5.625 x 10 3 , C f = 1.935 x 10~ 3 
Water flows between two parallel walls as shown in Fig. 9.14. The velocity is 
uniform at the entrance and core region. Beyond a distnace Le downstream from 
the entrance, the flow becomes fully developed so that the velocity varies over the 
entire width 2 h of the channel. In the boundary layer region, velocity varies as 

yV 

— where 8= a vx ; a being a constant. Determine the acceleration 
8 J 

on the axis of symmetry for 0 < x < Le. 


u = U(x) 
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T 

h 

1 


Fig. 9.14 Development of boundary layers in a channel 

Ans. a = (U 2 0 l Me) {xlLe)~ m l{ 1 - (2/3)(x/Z,e) 0 ' 5 ) 3 

9.9 Consider the laminar boundary layer on a flat plate with uniform suction velocity 
V 0 as shwon in Fig. 9.15: 



\ x I J J J J | 


Incompressible 

fluid 


= constant 


Fig. 9.15 Flow on a flat plate with uniform suction velocity 


Far down the plate (large x), a fully-developed situation may be shown to exist in 
which the velocity distribution does not vary with x. Find the velocity distribution 
in this region, as well as the wall shear. The governing equations are 


du 

dx 


dv 

dy 


= 0 


and 


du 

dx 


■ + v ■ 


dv 

dy 


dx: dy 


The boundary conditions are at y = 0, u = 0, v = V 0 and u(°°) = U 

Ans. u = U{ 1 -exp (-F 0 y/v)}, T„, = pUV 0 

9.10 Consider a laminar boundary layer on flat plate with a velocity profile given by 

— = — n - — n 3 , where, ri = v/8 
U 2 2 

For this profile — = ’ . 

x Re, 

Determine an expression for the local skin friction coefficient in terms of distance 
and flow properties. 

Ans. (C f = 0.647 (Re,) 0 ' 5 ) 

9.11 A low-speed wind-tunnel is provided with air supply upto a speed of 50 m/s at 
20 °C. One needs to study the behaviour of boundary-layer over a flat plate kept 
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inside the wind-tunnel, up to a Reynolds number of Re^ = 10 s . What is the 
minimum plate length that should be used? At what distance from the leading 
edge would the transition occur if the critical Reynolds number is Re,. c = 3.5 x 
10 5 ? At 25 °C, v of air is 15.7 x 1(T 6 m 2 /s. 


Ans. (x min = 31.4 m,x cr = 0.109 m) 


9.12 Perform a numerical solution (develop a FORTRAN code) using Eq. (9.24) and a 
Runge-Kutta method (as outlined in the text) which will iterate the Blasius 
equation from an initial guess H( 0) = 0.2 and converge to the exact value of 77(0) 
= 0.4696. 

9.13 A liquid film of uniform thickness flows down the inner wall of vertical pipe, the 
thickness of the film being very small compared to the pipe radius. Show that, in 
the absence of any appreciable shear force on the free surface of the film, the 
volume flow of liquid per unit time, Q x , is given by 


2 nrgt 2 


where, r is the pipe radius, g the gravitational acceleration, t is the thickness of 
the film and V is the kinematic viscosity of the liquid. 

Show further that, if air flows up the pipe at such a rate that the free surface 
of the film remains at rest, then the volume flow of liquid per unit time, Q, is given 
by 



4 v Pg dx 



dx 

and other symbols are as defined above. 

9.14 The velocity distribution in the laminar boundary layer is of the form 


u 


F(77)+AG(77) 


where, F(ri) = — t] - —; G{r\) = —- — + — and A = 

2 2 4 2 4 <5 


8 2 dUe 


v dx 


Find the value of A when the flow is on the point of separating and show that 
then the displacement thickness will be half the boundary layer thickness. 










The McGraw-Hill Companies 



Turbulent Flow 


10.1 INTRODUCTION 

The turbulent motion is an irregular motion. The irregularity associated with 
turbulence is such that it can be described by the laws of probability and turbulent 
fluid motion can be considered as an irregular condition of flow in which various 
quantities (such as velocity components and pressure) show a random variation 
with time and space in such a way that the statistical average of those quantities 
can be quantitatively expressed. 

An irregular motion is associated with random fluctuations. It is postulated 
that the fluctuations inherently come from disturbances (such as roughness of a 
solid surface) and they may be either dampened out due to viscous damping or 
may grow by drawing energy from the free stream. At a Reynolds number less 
than the critical, the kinetic energy of flow is not enough to sustain the random 
fluctuation against the viscous damping and in such cases laminar flow continues 
to exist. At somewhat higher Reynolds number than the critical Reynolds number, 
the kinetic energy of flow supports the growth of fluctuations and transition to 
turbulence takes place. 

10.2 CHARACTERISTICS OF TURBULENT FLOW 

In view of the preceding discussion, the most important characteristic of turbulent 
motion is the fact that velocity and pressure at a point fluctuate with time in a 
random manner (Fig. 10.1). 
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Up 


Up 



Laminar 


(a) 


Turbulent 


(b) 


Fig. 10.1 Variation of u components of velocity for laminar and 
turbulent flows at a point P 

The mixing in turbulent flow is more due to these fluctuations. As a result we 
can see more uniform velocity distributions in turbulent pipe flows as compared 
to laminar flows (see Fig. 10.2). 


Laminar Turbulent 




Fig. 10.2 Comparison of velocity profiles in a pipe for (a) laminar 
and (b) turbulent flows 


Turbulence can be generated by frictional forces at the confining solid walls or 
by the flow of layers of fluids with different velocities over one another. The 
turbulence generated in these two ways are considered to be different. Turbulence 
generated and continuously affected by fixed walls is designated as wall tur¬ 
bulence, and turbulence generated by two adjacent layers of fluid in absence of 
walls is termed as free turbulence. 

One of the effects of viscosity on turbulence is to make the flow more 
homogeneous and less dependent on direction. If the turbulence has the same 
structure quantitatively in all parts of the flow field, the turbulence is said to be 
homogeneous. Turbulence is called isotropic if its statistical features have no 
directional preference and perfect disorder persists. Its velocity fluctuations are 
independent of the axis of reference, i.e. invariant to axis rotation and reflection. 
Isotropic turbulence is by its definition always homogeneous. In such a situation, 
the gradient of the mean velocity does not exist. The mean velocity is either zero 
or constant throughout. However, when the mean velocity has a gradient, the 
turbulence is called anisotropic. 
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A little more discussion on homogeneous and isotropic turbulence is needed at 
this stage. The term homogeneous turbulence implies that the velocity fluctuations 
in the system are random. The average turbulent characteristics are independent 
of the position in the fluid, i.e., invariant to axis translation. 

Consider the root mean square velocity fluctuations 



In homogeneous turbulence, the rms values of u, V and w' can all be different, 
but each value must be constant over the entire turbulent field. It is also to be 
understood that even if the rms fluctuation of any component, say u s are constant 
over the entire field, the instantaneous values of u may differ from point to point 
at any instant. 

In addition to its homogeneous nature, if the velocity fluctuations are 
independent of the axis of reference, i.e., invariant to axis rotation and reflection, 
the situation leads to isotropic turbulence, which by definition as mentioned 
earlier, is always homogeneous. 

In isotropic turbulence fluctuations are independent of the direction of 
reference and 



or u — V = w' 

Again it is of relevance to say that even if the rms fluctuations at any point are 
same, their instantaneous values may differ from each other at any instant. 

Turbulent flow is also diffusive. In general, turbulence brings about better 
mixing of a fluid and produces an additional diffusive effect. The term “eddy- 
diffusion” is often used to distinguish this effect from molecular diffusion. The 
effects caused by mixing are as if the viscosity is increased by a factor of 100 or 
more. At a large Reynolds number there exists a continuous transport of energy 
from the free stream to the large eddies. Then, from the large eddies smaller eddies 
are continuously formed. Near the wall smallest eddies dissipate energy and 
destroy themselves. 

10.3 LAMINAR-TURBULENT TRANSITION 

For turbulent flow over a flat plate, the boundary layer starts out as laminar flow 
at the leading edge and subsequently, the flow turns into transition flow and very 
shortly thereafter turns into turbulent flow. The turbulent boundary layer 
continues to grow in thickness, with a small region below it called a viscous 
sublayer. In this sublayer, the flow is well behaved, just as the laminar boundary 
layer (Fig. 10.3). 

A careful observation further suggests that at a certain axial location, the 
laminar boundary layer tends to become unstable. Physically this means that the 
disturbances in the flow grow in amplitude at this location. 
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On 


U„ U QO 



Fig. 10.3 Laminar-turbulent transition 


Free stream turbulence, wall roughness and acoustic signals may be among the 
sources of such disturbances. Transition to turbulent flow is thus initiated with 
the instability in laminar flow. The possibility of instability in boundary layer 
was felt by Prandtl as early as 1912. The theoretical analysis of Tollmien and 
Schlichting showed that unstable waves could exist if the Reynolds number was 
575. The Reynolds number was defined as 

Re = U„tf/v 

where U,„ is the free stream velocity and <5 is the displacement thickness. Taylor 
developed an alternate theory, which assumed that the transition is caused by a 
momentary separation at the boundary layer associated with the free stream 
turbulence. In a pipe flow the initiation of turbulence is usually observed at 
Reynolds numbers ( U<JD/v ) in the range of 2000 to 2700. The development starts 
with a laminar profile, undergoes a transition, changes over to turbulent profile 
and then stays turbulent thereafter (Fig. 10.4). The length of development is of the 
order of 25 to 40 diameters of the pipe. The mechanisms related to the growth and 
the decay of turbulence in a flow field are indeed an advanced topic and beyond 
the scope of this text. However, the interested readers may like to refer to 
Tennekes and Lumley [1] and Hinze [2] for advanced knowledge. 

Developing 

^_ region " _*j_►pfully developed 

turbulent flow 



► .r 


Fig. 10.4 Development of turbulent flow in a circular duct 
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10.4 CORRELATION FUNCTIONS 

A statistical correlation can be applied to fluctuating velocity terms in turbulence. 
Turbulent motion is by definition eddying motion. Notwithstanding the 
circulation strength of the individual eddies, a high degree of correlation exists 
between the velocities at two points in space, if the distance between the points is 
smaller than the diameter of the eddy. 

Consider a statistical property of a random variable (velocity) at two points 
separated by a distance r. An Eulerian correlation tensor (nine terms) at the two 
points can be defined by 

Q = u(.v)u (x + r) 

In other words, the dependence between the two velocities at two points is 
measured by the correlations, i.e. the time averages of the products of the 
quantities measured at two points. The correlation of the u components of the 
turbulent velocity of these two points is defined as 

u'(x)u'(x + r ) 

It is conventional to work with the non-dimensional form of the correlation, such 
as 

= u'{x)u'{x + r) 

— — 1/2 

\u' 2 ( x)u' 2 ( x 4- r)J 

A value of R(r ) of unity signifies a perfect correlation of the two quantities 
involved and their motion is in phase. Negative value of the correlation function 
implies that the time averages of the velocities in the two correlated points have 
different signs. Figure 10.5 shows typical variations of the correlation R with 
increasing separation r. 

To describe the evolution of a fluctuating function u'(t), we need to know the 
manner in which the value of u at different times are related. For this purpose a 
correlation function 

R(z) = u'{t)u'(t + t)/u' 2 

between the values of u' at different times has been chosen. This is called 
autocorrelation function. 

The correlation studies reveal that the turbulent motion is composed of eddies 
which are convected by the mean motion. The eddies have a wide range variation 
in their size. The size of the large eddies is comparable with the dimensions of the 
neighbouring objects or the dimensions of the flow passage. The size of the 
smallest eddies can be of the order of 1 mm or less. However, the smallest eddies 
are much larger than the molecular mean free paths and the turbulent motion 
obeys the principles of continuum mechanics. 
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R (r) k 1.0 


Fig. 10.5 Variation of R with the distance of separation r 


10.5 MEAN MOTION AND FLUCTUATIONS 


In 1883, O. Reynolds conducted experiments with pipe flow by feeding into the 
stream a thin thread of liquid dye. For low Reynolds numbers, the dye traced a 
straight line and did not disperse. With increasing velocity, the dye thread got 
mixed in all directions and the flowing fluid appeared to be uniformly coloured in 
the downstream. It was conjectured that on the main motion in the direction of the 
axis, there existed a superimposed motion all along the main motion at right 
angles to it. The superimposed motion causes exchange of momentum in 
transverse direction and the velocity distribution over the cross-section is more 
uniform than in laminar flow. This description of turbulent flow which consists of 
superimposed streaming and fluctuating (eddying) motion is referred to as 
Reynolds decomposition of turbulent flow. 

Here, we shall discuss different descriptions of mean motion. Generally, for 
Eulerian velocity u, the following two methods of averaging could be obtained, 
(i) Time average for a stationary turbulence 



(ii) Space average for a homogeneous turbulence 



-X 


For a stationary and homogeneous turbulence, it is assumed that the two 
averages lead to the same result: u* = u s and the assumption is known as the 
ergodic hypothesis. 

In our analysis, average of any quantity will be evaluated as a time average. 
We take t x a finite interval. This interval must be larger than the time scale of 
turbulence. Needless to say that it must be small compared with the period t 2 of 
any slow variation (such as periodicity of the mean flow) in the flow field that we 
do not consider to be chaotic or turbulent. 

Thus, for a parallel flow, it can be written that the axial velocity component is 


u(y, t) = u(y) + u (r, 0 


(10.1) 











The McGraw-Hill Companies 


404 | |Introduction to Fluid Mechanics and Fluid Machines 


As such, the time mean component u (y) determines whether the turbulent 
motion is steady or not. Let us look at two different turbulent motions described in 
Fig. 10.6 (a) and (b). The symbol T signifies any of the space variables. 

While the motion described by Fig. 10.6 (a) is for a turbulent flow with steady 
mean velocity the Fig. 10.6 (b) shows an example of turbulent flow with unsteady 
mean velocity. The time period of the high frequency fluctuating component is t l 
whereas the time period for the unsteady mean motion is t 2 and for obvious reason 
t 2 »t x . Even if the bulk motion is parallel, the fluctuation u being random varies 
in all directions. Now let us look at the continuity equation, 

du dv dw 

-+-+-= 0 

dx dy dz 

Invoking Eq. (10.1) in the above expression, we get 


du du' dv dw 

- 1 - 1 - 1 - 

dx dx dy dz 


( 10 . 2 ) 



(a) Steady (b) Unsteady 

Fig. 10.6 Steady and unsteady mean motions in a turbulent flow 


Since ^ 0, Eq. (10.2) depicts that y and z components of velocity exist even 
dx 

for the parallel flow if the flow is turbulent. We can write 


u(y, t) = u (y) + u (r, t) 
v = o + v' (r, t) 
w — o + w' (r, t ) 


(10.3) 


However, the fluctuating components do not bring about the bulk displacement 
of a fluid element. The instantaneous displacement is u d t, and that is not 
responsible for the bulk motion. We can conclude from the above 


T 

| u d t = 0, (t| < T < t 2 ) 

-T 

Due to the interaction of fluctuating components, macroscopic momentum 
transport takes place. Therefore, interaction effect between two fluctuating 
components over a long period is non-zero and this can be expressed as 
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J u v' dt ^ 0 

-T 

We take time average of these two integrals and write 


, J 

u' = — [ u dt = 0 
IT J 


and 


2 T 

— f u v' dt± 0 

IT J 


-T 

T 


IT 


(10.4a) 


(10.4b) 


-T 


Now, we can make a general statement with any two fluctuating parameters, say, 
with/' and g as 

T=r= 0 (10.5a) 

TV^O (10.5b) 

The time averages of the spatial gradients of the fluctuating components also 
follow the same laws, and they can be written as 


df _ d 2 r _ 


ds 


ds 2 


= 0 


and 


d(f'g') 


ds 


*0 


( 10 . 6 ) 


The intensity of turbulence or degree of turbulence in a flow is described by 
the relative magnitude of the root mean square value of the fluctuating 
components with respect to the time averaged main velocity. The mathematical 
expression is given by 

I=Jj(^ 2 +^ 2 +^ 2 ) lu„ (10.7a) 

The degree of turbulence in a wind tunnel can be brought down by introducing 
screens of fine mesh at the bell mouth entry. In general, at a certain distance from 
the screens, the turbulence in a wind tunnel becomes isotropic, i.e. the mean 
oscillation in the three components are equal, 


In this case, it is sufficient to consider the oscillation u in the direction of flow 
and to put 


7=V^"/l7 00 (10.7b) 

This simpler definition of turbulence intensity is often used in practice even in 
cases when turbulence is not isotropic. 

Following Reynolds decomposition, it is suggested to separate the motion into 
a mean motion and a fluctuating or eddying motion. Denoting the time average of 
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the u component of velocity by u and fluctuating component as u, we can write 
down the following, 

u = u + u, v = V + v', w = w + w', p = p + p ( 10 . 8 ) 

By definition, the time averages of all quantities describing fluctuations are 
equal to zero. 

* 7 = 0 , 7 = 0 , >7 = 0 , 7 = 0 ( 10 . 9 ) 

The fluctuations u, v', and w' influence the mean motion it , v and w in such 
a way that the mean motion exhibits an apparent increase in the resistance to 
deformation. In other words, the effect of fluctuations is an apparent increase in 
viscosity or macroscopic momentum diffusivity. 

We shall state some rules of mean time-averages here. If / and g are two 
dependent variables and if s denotes any one of the independent variables x, y, z, 
t, then 

/ = /; f + g = f + g',f g = f g’, 

J 7 ^ = j 7 ds (lo.io) 


10.6 DERIVATION OF GOVERNING 

EQUATIONS FOR TURBULENT FLOW 

For incompressible flows, the Navier-Stokes equations can be rearranged in the 
form 


du d(u 2 ) d(uv) d(uw ) 
dt dx dy dz 

(9v d(uv ) d(v 2 ) d(vw) 

dt dx dy dz 

dw d(uw) d(vw) <9(w 2 ) 
dt dx dy dz 


= - + p v 2 u 

dx 


dp 

dy 

d£ 

dz 


pV 2 \ 


■ /uV 2 w 


and 


w 


(10.11a) 

(10.11b) 

(10.11c) 

( 10 . 12 ) 


du dv „ „ , 

-+-+-= 0 

dx dy dz 

Let us express the velocity components and pressure in terms of time-mean 
values and corresponding fluctuations. In continuity equation, this substitution 
and subsequent time averaging will lead to 


d{u + u') + d(v 4- v') + d(w + w') _ 


or 


dx dy 

^du dv dw 

- 1 - 1 - 

dx dy dz 


dz 

f du' dv' dw 


dx 


• + • 
dy dz 


= 0 
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Since = dv = dw = Q [From Eq. (10.6)] 

dx ay dz 

We can write 


du dv dw 

- 1 - 1 - 

dx dy dz 

From Eqs (10.13a) and (10.12), we obtain 


(10.13a) 


du' dv' dw' 

- 1 - 1 - 

dx dy dz 


(10.13b) 


It is evident that the time-averaged velocity components and the fluctuating 
velocity components, each satisfy the continuity equation for incompressible 
flow. Let us imagine a two-dimensional flow in which the turbulent components 
are independent of the z-direction. Eventually, Eq. (10.13b) tends to 


du' _ dv' 
dx dy 


(10.14) 


On the basis of condition (10.14), it is postulated that if at an instant there is an 
increase in u in the x-direction, it will be followed by an increase in v' in the 

negative y-direction. In other words, u'v' is non-zero and negative. This is 
another important consideration within the framework of mean-motion 
description of turbulent flows. 

Invoking the concepts of (10.8) into the equations of motion (10.1 la, b, c), we 
obtain expressions in terms of mean and fluctuating components. Now, forming 
time averages and considering the rules of (10.10) we discern the following. The 


terms which are linear, such as 


du' 

~dT 


and 


d 2 u' 

dx 2 


vanish when they are averaged 


[from (10.6)]. The same is true for the mixed terms like u • u , or u • v', but the 
quadratic terms in the fluctuating components remain in the equations. After 

averaging, they form u ' 2 , u'v' etc. 

For example, if we perform the aforesaid exercise on the x momentum 
equation, we shall obtain 


P 


d(u + u') d(u + u'Y d(u + w')(v + v') 
dt + + 


d(p + p ') 
dx 


dx 

dy 

+ P 

d 2 (u +u') d 

dx 2 


2 /— 


dy 


d(u + u')(w+ w') 
dz 


d 2 (u +u') 

+ —^—- 
dz 2 


or 



d(u 1 +u' z ) d(u ■ v + u'v') d(u ■ w + u'w') 
+ -+-^ + - 


dx 


dy 


dz 
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~)2 — ~\2 — -) 2 — ( ~\2—j 

a u o u o u a 


; o "*- T 

dx 2 dy 2 dz 2 


d/X 


° 

o uf 


< 2-2 \ 
O UT 


or 


d(u") d(u-v) d(u -w) 


dx 


dy 


dz 


— — ~~~ + n V” u 
dx 


d — d 

-M V +- 


dx dy dz 

Introducing simplifications arising out of continuity Eq. (10.13a), we shall obtain 

dp 


_ du _ du _ du 

U -h V-b W - 

dx dy dz 


~P 


dx 


dy 


dx 


dz 


- /J. V - u 


Performing a similar treatment on y and z momentum equations, finally we 
obtain the momentum equations in the form 


_ du _ du _ du 

U -b V-b W - 

dx dy dz 


- pV 2 it 


~P 


.'2 


dx 


dy 


_<9v _<9v _<9v 

U -b V-b W - 

dx dy dz 


dp 

dx 

d_ 

dz 


dy 


(10.15a) 


d —r~. dp ^ . 

-u v + — v + — v w 


dx 


dy 


P 

_ dw _ dw _ dw 

U -b V-b W - 

dx dy dz 


d_ 

dz 


(10.15b) 


dp 

dz 


+ pV 2 w 


dx 


dy 


dz 


-w 


(10.15c) 


The left hand side of Eqs (10.15a)-(10.15c) are essentially similar to the 
steady-state Navier-Stokes equations if the velocity components u, v and w are 
replaced by u , V and w. The same argument holds good for the first two terms on 
the right hand side of Eqs (10.15a)-(10.15c). However, the equations contain 
some additional terms which depend on turbulent fluctuations of the stream. These 
additional terms can be interpreted as components of a stress tensor. Now, the 
resultant surface force per unit area due to these terms may be considered as 


_ du _ du _ du 

U — -b v—-b — 

dx dy dz 


— - + pV 2 u 

dx 
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d , d , d , 

-<7 -\ -T H-T 

xx -j yx -i zx 

dx dy dz 

= - — + ,uV 2 ro 

dy 

a , 

7 4-T 

^ <?z ^ 


(10.16a) 


_ <?G3 — <905 _ <905 

u -1- 05-1- w - 

dx dy dz 


i- T ' 


<?v 


(10.16b) 


_ dw _ <9w _ dw 

U — -h v—-h Vl ; —— 

dx dy dz 


=- +iuv 2 w 

dz 


+ 


d . , d d 

Ttz -) -j 

dx dy dz 


(10.16c) 


Comparing Eqs (10.15) and (10.16), we can write 


1 

T' 

L xy 

_ 1 


/2 

U 

U V 

1 

> 

<y 

°yy 


= -P 

11 V 

V' 2 

v'w' 

<z 




u'w' 

v'vri 

/2 

w 


(10.17) 


It can be said that the mean velocity components of turbulent flow satisfy the 
same Navier-Stokes equations of laminar flow. However, for the turbulent flow, 
the laminar stresses must be increased by additional stresses which are given by 
the stress tensor (10.17). These additional stresses are known as apparent stresses 
of turbulent flow or Reynolds stresses. Since turbulence is considered as eddying 
motion and the aforesaid additional stresses are added to the viscous stresses due 
to mean motion in order to explain the complete stress field, it is often said that 
the apparent stresses are caused by eddy viscosity. The total stresses are now 


- du , 2 

= - p + 2 p — - pu 

dx 


f du 


= M 


<?v A 


+ 

dy dx 


- pu v 


(10.18) 


and so on. The apparent stresses are much larger than the viscous components, 
and the viscous stresses can even be dropped in many actual calculations. 


10.7 TURBULENT BOUNDARY LAYER 
EQUATIONS 


For a two-dimensional flow (w = 0) over a flat plate, the thickness of turbulent 
boundary layer is assumed to be much smaller than the axial length and the order 
of magnitude analysis (refer to Chapter-9) may be applied. As a consequence, the 
following inferences are drawn: 


(a) 


dp_ = 
dy 


0 , 



d p_ 

dx 
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(c) 


d 2 ~ 


dx 


u 

T « 


d 2 - 


dy 


u , 
, , and 


<d) ic(- pu ' 2 ) << i r (- puV '> 


dy 

The turbulent boundary layer equation together with the equation of continuity 
becomes 


du dv . 

-+-= 0 

dx dy 

_ dii - du 1 d p d 

U -b V --- V —— 

dx dy p dx dy 


du —f—, 
v— - U V 


dy 


(10.19) 


( 10 . 20 ) 


A comparison of Eq. (10.20) with laminar boundary layer Eq. (9.10) depicts, 
that: u, v and p are replaced by the time average values u , v and p , and laminar 

viscous force per unit volume is replaced by (r t + T f ), where x x =n^- 

dy dy dy 

is the laminar shear stress and T t = - pit' v' is the turbulent stress. 

10.8 BOUNDARY CONDITIONS 


All the components of apparent stresses vanish at the solid walls and only stresses 
which act near the wall are the viscous stresses of laminar flow. The boundary 
conditions, to be satisfied by the mean velocity components, are similar to laminar 
flow. A very thin layer next to the wall behaves like a near wall region of the 
laminar flow. This layer is known as laminar sublayer and its velocities are such 
that the viscous forces dominate over the inertia forces. No turbulence exists in it 
(see Fig. 10.7). For a developed turbulent flow over a flat plate, in the near wall 
region, inertial effects are insignificant, and we can write from Eq. (10.20), 


d 2 U d{u' V') _ 

dy 2 dy 



Fig. 10.7 Different zones of a turbulent flow past a wall 
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which can be integrated as,- u' v' = constant 

dy 

Again, as we know that the fluctuating components, do not exist near the wall, 
the shear stress on the wall is purely viscous and it follows 


v 


du 

dy 


y= o 


T 


W 


P 


However, the wall shear stress in the vicinity of the laminar sublayer is estimated 
as 


Tv = P 


Us~ o 
<5-0 


U. 


(10.21a) 


where U s is the fluid velocity at the edge of the sublayer. The flow in the sublayer 
is specified by a velocity scale (characteristic of this region). We define the 
friction velocity, 


1/2 


(10.21b) 


as our velocity scale. Once u T is specified, the structure of the sublayer is 
specified. It has been confirmed experimentally that the turbulent intensity 
distributions are scaled with u t . For example, maximum value of the u' 2 is 
always about 8 u 2 . The relationship between w T and the U s can be determined from 
Eqs (10.21a) and (10.21b) as 

2 U s 

u,= v— 1L 

T <5, 

Let us assume (7 = CT/. Now we can write 


2 - U„ - 

u~ = Cv —— where C is a proportionality constant (10.22a) 


or 


S s»T = £ 


(10.22b) 


y n 

Hence, a non-dimensional coordinate may be defined as, 1 ) = -—— which will 

v 

help us estimating different zones in a turbulent flow. The thickness of laminar 
sublayer or viscous sublayer is considered to be ij ~ 5. Turbulent effect starts in 
the zone of t] > 5 and in a zone of 5 < T) < 70, laminar and turbulent motions 
coexist. This domain is termed as buffer zone. Turbulent effects far outweigh the 
laminar effect in the zone beyond rj = 70 and this regime is termed as turbulent 
core. 

For flow over a flat plate, the turbulent sheat stress {-pu'v') is constant 
throughout in the y direction and this becomes equal to T w at the wall. In the event 
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of flow through a channel, the turbulent shear stress (- pu'v') varies with y and 
it is possible to write 


= 


T 


w 


I 

h 


(10.22c) 


where the channel is assumed to have a height 2 h and C is the distance measured 
from the centreline of the channel (= h - y ). Figure 10.7 explains such variation of 
turbulent stress. 


10.9 SHEAR STRESS MODELS 


In analogy with the coefficient of viscosity for laminar flow, J. Boussinesq 
introduced a mixing coefficient p z for the Reynolds stress term by invoking 


— r~, du 

T t = - pu v =p t — 
dy 

Now the expressions for shearing stresses are written as 


du du 

T /=PV-^,T,=7b — 
dy dy 


pv 


‘ dy 


such that the equation 


_ du - du 

U lh +V lfy 

may be written as 


_ du - du 

U lh +V lfy 


1 dp d du ; 
~p'^c + ~dy[~d^~ U V 


1 dp d 
p dx dy 


j(v + v,) 



(10.23) 


The term v, is known as eddy viscosity and the model is known as eddy 
viscosity model. The difficulty in using Eq. (10.23) can be discussed herein. The 
value of V, is not known. The term v is a property of the fluid whereas V, is 
attributed to random fluctuations and is not a property of the fluid. However, it is 
necessary to find out empirical relations between v t and the mean velocity. We 
shall discuss one such well known relation between the aforesaid apparent or 
eddy viscosity and the mean velocity components in the following subsection. 


10.9.1 Prandtl's Mixing Length Hypothesis 

Let us consider a fully developed turbulent boundary layer (Fig. 10.3). The 
streamwise mean velocity varies only from streamline to streamline. The main 
flow direction is assumed parallel to the x-axis (Fig. 10.8). 

The time average components of velocity are given by u = u (y), V = 0, 
w = 0. The fluctuating component of transverse velocity v' transports mass and 
momentum across a plane at y 1 from the wall. The shear stress due to the 
fluctuation is 
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given by 


T = 


—r~, du 
pu v =P,— 
dy 


(10.24) 


A lump of fluid, which comes to the layerjh from a layer (y x — /) has a positive 
value of v'. If the lump of fluid retains its original momentum then its velocity at 
its current location^ is smaller than the velocity prevailing there. The difference 
in velocities is then 


A u | 


it Oh) - u Oh - 0 = l 


f du^ 


dy 


(10.25) 



Fig. 10.8 One-dimensional parallel flow and Prandtl's mixing length hypothesis 


The above expression is obtained by expanding the function u Oh - l) in a 
Taylor series and neglecting all higher order terms and higher order derivatives. 
As it is said, / is a small length scale known as Prandtl’s mixing length. Prandtl 
proposed that the transverse displacement of any fluid particle is, on an average, 
7’. Let us consider another lump of fluid with a negative value of v'. This is 
arriving aty 1 from (y \ + /)• If this lump retains its original momentum, its mean 
velocity at the current lamina y 1 will be somewhat more than the original mean 
velocity ofV[. This difference is given by 


Am 2 


u (jh + /) _ « (Jh) ~ 1 


f du} 

\ d yj 


(10.26) 


The velocity differences caused by the transverse motion can be regarded as the 
turbulent velocity components at y 1 . We calculate the time average of the absolute 
value of this fluctuation as 


\ u '\ - ~ (| Ahj | +1 Aw 2 |) = / 


( du^ 


\ d yj 


(10.27) 
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Suppose these two lumps of fluid meet at a Iaycry,. The lumps will collide with a 
velocity 2 u and diverge. This proposes the possible existence of transverse 
velocity component in both directions with respect to the layer at y x . Now, 
suppose that the two lumps move away in a reverse order from the 1 ayer y, with a 
velocity 2 u . The empty space will be filled from the surrounding fluid creating 
transverse velocity components which will again collide at y ,. Keeping in mind 
this argument and the physical explanation accompanying Eqs (10.14), we may 
state that 


\ v '\ 


or \v’\ = const \u'\ *= (const) / 

along with the condition that the moment at which u is positive, v' is more likely 
to be negative and conversely when u is negative. Possibly, we can write at this 
stage 

u'v' = - Cj \u f \ \v'\ 


du 

dy 


or 


u v 


c 2 i l 


r du^ 

, d y. 


(10.28) 


where Cj and C 2 are different proportionality constants. However, the constant 
C 2 can now be included in still unknown mixing length and Eq. (10.28) may be 
rewritten as 


u'v' =-l 2 


du 

\ d yj 


For the expression of turbulent shearing stress T r , we may write 

du 


T t = -p u'v' = pi 


\dy j 


(10.29) 


After comparing this expression with the eddy viscosity concept and 
Eq. (10.24), we may arrive at a more precise definition, 


T, = pl~ 




du 

‘~dy 

where the apparent viscosity may be expressed as 


du 

fdu' 

dy 

UtJ 


MrPf 


du 


dy 


and the apparent kinematic viscosity is given by 


v t = l 2 


du 

dy 


(10.30a) 


(10.30b) 


(10.30c) 
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The decision of expressing one of the velocity gradients ofEq. (10.29) in terms 
du 


of its modulus as 


dy 


was made in order to assign a sign to T, according to the 


sign of-. It may be mentioned that the apparent viscosity and consequently, 

dy 

the mixing length are not properties of fluid. They are dependent on turbulent 
fluctuation. However, our problem is still not resolved. How to determine the 
value of the mixing length? Several correlations, using experimental results 
for T t have been proposed to determine l. 

However, so far the most widely used value of mixing length in the regime of 
isotropic turbulence is given by 

l = X y (10.31) 

where y is the distance from the wall and X is known as von Karman constant 
0 = 0.4). 


10.10 UNIVERSAL VELOCITY DISTRIBUTION LAW 

AND FRICTION FACTOR IN DUCT FLOWS FOR 
VERY LARGE REYNOLDS NUMBERS 


For flows in a rectangular channel at very large Reynolds numbers the laminar 
sublayer can practically be ignored. The channel may be assumed to have a width 
2 h and the x axis will be placed along the bottom wall of the channel. We shall 
consider a turbulent stream along a smooth flat wall in such a duct and denote the 
distance from the bottom wall by y, while u(y ) will signify the velocity. In the 
neighbourhood of the wall, we shall apply 

i = xy 


According to Prandtl’s assumption, the turbulent shearing stress will be 


t = p/y 2 


( du} 

\ d yj 


(10.32) 


At this point, Prandtl introduced an additional assumption which like a plane 
Couette flow takes a constant shearing stress throughout, i.e 

t, = t w (10.33) 

where T w denotes the shearing stress at the wall. Invoking once more the friction 

-|l/2 


velocity w = 


, we obtain 


uhx 1 / 


f 3-\ 2 

du 


\ d yj 


(10.34) 


or 


du _ u T 
dy Xy 


On integrating we find 


(10.35) 
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II 

u = —— In y + C (10.36) 

X 

Despite the fact that Eq. (10.36) is derived on the basis of the friction velocity in 
the neighbourhood of the wall because of the assumption that T w = T, = constant, 
we shall use it for the entire region. At y = h (at the horizontal mid plane of the 
channel), we have u = [/ max . The constant of integration is eliminated by 
considering 


U max = ^lnh + C 

X 


or 


C = U max — In h 


T 

X 

Substituting C in Eq. (10.36), we get 

f ,\ 


and U 'max « = 1 ln 

U T X 


\y) 


(10.37) 


Equation (10.37) is known as universal velocity defect law of Prandtl and its 
distribution has been shown in Fig. 10.9. 

Flere, we have seen that the friction velocity u t is a reference parameter for 
velocity. We shall now discuss the problem with (u lu t ) and rj (= y u T /v ) as 
parameters. Equation (10.36) can be rewritten once again for this purpose as 

— = — ln v + C 

«T X 



0 0.4 0.8 1.0 

yl h 


Fig. 10.9 Distribution of universal velocity defect law of 
Prandtl in a turbulent channel flow 

The no-slip condition at the wall cannot be satisfied with a finite constant of 
integration. This is expected that the appropriate condition for the present problem 
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should be that u = 0 at a very small distance y = y (t from the wall. Hence, 
Eq. (10.36) becomes 

— = ' (1 n y - In y 0 ) (10.38) 

11 T X 

The distance y 0 is of the order of magnitude of the thickness of the viscous layer. 
Now we can write Eq. (10.38) as 

— = — In y — - In /3 

u t X L V 

or — =A x \mj + D x (10.39) 

U T 

where A x = (\/X), the unknown /j is included in D x . 

Equation (10.39) is generally known as the universal velocity profile because 
of the fact that it is applicable from moderate to a very large Reynolds number. 
However, the constants A x and D x have to be found out from experiments. The 
aforesaid profile is not only valid for channel (rectangular) flows, it retains the 
same functional relationship for circular pipes as well. It may be mentioned that 
even without the assumption of having a constant shear stress throughout, the 
universal velocity profile can be derived. Interested readers are referred to 
Example 10.3. 

Experiments, performed by J. Nikuradse, showed that Eq. (10.39) is in good 
agreement with experimental results. Based on Nikuradse’s and Reichardt’s 
experimental data, the empirical constants of Eq. (10.39) can be determined/or a 
smooth pipe as 

— = 2.5 In 7] + 5.5 (10.40) 

u T 

This velocity distribution has been shown through curve (b) in Fig. 10.10. 



sublayer zone zone 


Fig. 10.10 The universal velocity distribution law for smooth pipes 
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However, the corresponding friction factor concerning Eq. (10.40) is 

= 2.0 log 10 (Re if ) - 0.8 (10.41) 

As mentioned earlier, the universal velocity profile does not match very close to the 
wall where the viscous shear predominates the flow. However, von Karman 
suggested a modification for the laminar sublayer and the buffer zone which are 


if 


— = rj= —— for rj < 5.0 


(10.42) 


— = 11.5 log 10 ^ -3.0 for 5 < 77 < 60 (10.43) 

u t v 

Equation (10.42) has been shown through curve (a) in Fig. 10.10. 

It may be worthwhile to mention here that a surface is said to be hydraulically 
smooth so long 

£ „ u T 

0 < —2—L < 5 (10.44) 

v 

where e p is the average height of the protrusions inside the pipe. 

Physically, the above expression means that for smooth pipes protrusions will 
not be extended outside the laminar sublayer. If protrusions exceed the thickness of 
laminar sublayer, it is conjectured (also justified though experimental verification) 
that some additional frictional resistance will contribute to pipe friction due to the 
form drag experienced by the protrusions in the boundary layer. In rough pipes 
experiments indicate that the velocity profile may be expressed as: 

— = 2.5 In + 8.5 (10.45) 

U T £p 

At the centre-line, the maximum velocity is expressed as 

-^L = 2.5 In — + 8.5 (10.46) 

U z £ P 

Note that v no longer appears with R and £ p . This means that for completely rough 
zone of turbulent flow, the profile is independent of Reynolds number and a 
strong function of pipe roughness. However, for pipe roughness of varying 
degrees, the recommendation due to Colebrook and White works well. Their 
formula is 


1 

if 


1.74-2.0 log 10 



18.7 

ReV7 


(10.47) 


where R is the pipe radius. 

For £ p —> 0, this equation produces the result of the smooth pipes (Eq. (10.41)). 
For Re —» it gives the expression for friction factor for a completely rough pipe 
at a very high Reynolds number which is given by 
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on 


/= 7 --- -j (10.48) 

( R X 

2 log—+ 1.74 

v £ P J 

Turbulent flow through pipes has been investigated by many researchers because 
of its enormous practical importance. 

In the next section, we shall discuss, in detail the velocity distribution and 
other important aspects of turbulent pipe flows. 

10.11 FULLY DEVELOPED TURBULENT FLOW IN A 
PIPE FOR MODERATE REYNOLDS NUMBERS 


The entry length of a turbulent flow is much shorter than that of a laminar flow, J. 
Nikuradse determined that a fully developed profile for turbulent flow can be 
observed after an entry length of 25 to 40 diameters. We shall focus herein our 
attention to fully developed turbulent flow. Considering a fully developed 
turbulent pipe flow (Fig. 10.11) we can write 

2 71R z w = - j^J kR 1 (10.49) 


or 


dp) = 
dx J R 


(10.50) 



Fig. 10.11 Fully developed turbulent pipe flow 


It can be said that in a fully developed flow, the pressure gradient balances the 
wall shear stress only and has a constant value at any x. However, the friction 
factor (Darcy friction factor) is defined in a fully developed flow as 


PfUl 

2D 


(10.51) 


Comparing Eq. (10.50) with Eq. (10.51), we can write 

\=jPUl (10.52) 

H. Blasius conducted a critical survey of available experimental results and 
established the empirical correlation for the above equation as 
/= 0.3164 Re~° 25 , where Re = pU^D/p 


(10.53) 
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It is found that the Blasius ’ s formula is valid in the range of Reynolds number of 
Re < 10 5 . At the time when Blasius compiled the experimental data, results for 
higher Reynolds numbers were not available. However, later on, J. Nikuradse 
carried out experiments with the laws of friction in a very wide range of Reynolds 
numbers, 4x 10 3 < Re< 3.2 x 10 (> . The velocity profile in this range follows: 


u 

il 


1 In 


(10.54) 


where u is the time mean velocity at the pipe centre an dp is the distance from the 
wall. The exponent n varies slightly with Reynolds number. In the range of 
Re ~ 10 5 , n is 7. 

The ratio of u and U in for the aforesaid profile is found out by considering the 
volume flow rate Q as 

R 

Q = kR 2 U sw =^ 2 k r u dr 


or 


kR 2 U, v = 2ku J (R-y)(y/R) Vn (-dy) 


or 


or 


or 


or 


nR~ U„, = 2n i 


kR 2 U„, = 2k u 


( 


n 


-1 n +1 


n +1 
2 n 


R n y ” 

V J 


2/7+1 


( 2n + l _1> 

y n R " 


V 


J 


R 


H R 2 


kR 2 U„=2kR 2 F 


U.„ 


n +1 2«+l 

n 2 

(n +1) (2/7 + 1) 


(10.55a) 


u (n +1) (2/7 + 1) 

Now, for different values of n (for different Reynolds numbers) we shall obtain 
different values of U m /u from Eq. (10.55a). On substitution of Blasius resistance 
formula (10.53) in Eq. (10.52), the following expression for the shear stress at the 
wall can be obtained. 

0.3164 


T, = 


r> ^0.25 _ T r2 

Re p U m 


/ 


\ 1/4 


or 


or 


or 


T w = 0.03955 pU~ 


2RU... 


V““ “av / 
\l/4 


Tv= 0.03325 pUl' 4 (J 

n7/4 


t w = 0.03325 p 


U. 


av 

u y 


(m) 


7/4 ( y _ 

R 


\l/4 
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For n = 7, U, lv /u becomes equal to 0.8. Substituting U w !u = 0.8 in the above 
equation, we get 


T w = 0.03325 p (0.8 ) 1/4 (m ) 7/4 (v/R) m 

Finally it produces T w = 0.0225 p (u) 7/4 (v/R) 114 (10.55b) 


or 


u A x = 0.0225 (u) m | - 


1/4 


where u T is friction velocity. Flowever, u x may be splitted into « 7/4 and ii' z 4 and 
we obtain 


( - \ 
u 

7/4 , . 

( u t R\ 


= 44.44 


\ v ) 


or 


— = 8.74 
u T 


u r R 


(10.56a) 


Now we can assume that the above equation is not only valid at the pipe axis 
(y = R) but also at any distance from the wal 1 y and a general form is proposed as 

/ \l/7 

— = 8.74 (10.56b) 

« T V v ) 


In conclusion, it can be said that (l/7)th power velocity distribution law (10.56b) 
can be derived from Blasius ’ s resistance formula (10.53). Equation (10.55b) gives 
the shear stress relationship in pipe flow at a moderate Reynolds number, i.e Re< 
10 5 . Unlike very high Reynolds number flow, here laminar effect cannot be 
neglected and the laminar sublayer brings about remarkable influence on the outer 
zones. 

It is worth mentioning that the friction factor for pipe flows, /, defined by 
Eq. (10.53) is valid for a specific range of Reynolds number and for a particular 
surface condition. The experimental results for a wide range of Reynolds numbers 
and variety of pipe roughness can be summarized through Moody diagram which 
has been shown in Chapter 11. 


10.12 SKIN FRICTION COEFFICIENT FOR 

BOUNDARY LAYERS ON A FLAT PLATE 


Calculations of skin friction drag on lifting surface and on aerodynamic bodies are 
somewhat similar to the analyses of skin friction on a flat plate. Because of zero 
pressure gradient, the flat plate at zero incidence is easy to consider. In some of the 
applications cited above, the pressure gradient will differ from zero but the skin 
friction will not be dramatically different so long there is no separation. 

We begin with the momentum integral equation for flat plate boundary layer 
which is valid for both laminar and turbulent flow. 


_d_ 

dx 



T 


w 


p 


(10.57a) 
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Invoking the definition of C 


fa 


C r = 


ft ~ 1 2 

0 P U ~ , 
2. y 


C fa = 2 


dd 

dx 


, Eq. (10.57a) can be rewritten as 


(10.57b) 


Due to the similarity in the laws of wall, correlations of previous section may be 
applied to the flat plate by substituting 8 fori? and U m for the time mean velocity at 
the pipe centre. The rationale for using the turbulent pipe flow results in the situation 
of a turbulent flow over a flat plate is to consider that the time mean velocity, at the 
centre of the pipe is analogous to the free stream velocity, both the velocities being 
defined at the edge of boundary layer thickness. 

Finally, the velocity profile will be [following Eq. (10.54)] 


u 

U 


1/7 


for Re < 10' 


If we evaluate momentum thickness with this profile, we shall obtain 




<5/ \l/7 

y 


i-l4 


fan 


dv = —<5 
' 72 


(10.58) 


(10.59) 


Consequently, the law of shear stress (in range of Re < 10 5 ) for the flat plate is 
found out by making use of the pipe flow expression of Eq. (10.55b) as 


T w = 0.0225 


\ 1/4 


or 


—\2 


p(u) 


= 0.0225 


Ru 


1/4 


Substituting for u and 8 for R in the above expression, we get 

-il/4 


or 


put 


= 0.0225 


< 5£7 


(10.60) 


Once again substituting Eqs (10.59) and (10.60) in Eq. (10.57), we obtain 

_ r -]l/4 

— •— = 0.0225 —— 

72 dx- 8U 


or 


or 


Si/ 4 — = 02314 

dx 

<5 S/4 = 0.2892 


1/4 


u 


f v V/4 


\ U -J 


c 


(10.61) 
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For simplicity, if we assume that the turbulent boundary layer grows from the 
leading edge of the plate we shall be able to apply the boundary conditions x = 0, 
5=0 which will yield C = 0, and Eq. (10.61) will become 

r -ii/4 


(<S/.r) 5/4 = 0.2892 


v 

xU„ 


or 


or 

where 


S_ 

x 


0.37 



— = 0.37 (Re r r 1/5 
x 

Re v = (U r „ x)/v 


(10.62) 


From Eqs (10.57b), (10.59) and (10.62), it is possible to calculate the average 
skin friction coefficient on a flat plate as 

C f = 0.072 (Re z r 1/5 (10.63) 

It can be shown that Eq. (10.63) predicts the average skin friction coefficient 
correctly in the regime of Reynolds number below 2x 10 6 . 

This result is found to be in good agreement with the experimental results in the 
range of Reynolds number between 5 x 10 5 and 10 7 which is given by 

Cf = 0.074 (R e/ )r 1/5 (10.64) 

Equation (10.64) is a widely accepted correlation for the average value of 
turbulent skin friction coefficient on a flat plate. 

With the help of Nikuradse’s experiments, Schlichting obtained the semi- 
empirical equation for the average skin friction coefficient as 


C / = 


0.455 

(log Re) 2 ' 58 


(10.65) 


Equation (10.65) was derived assuming the flat plate to be completely turbulent 
over its entire length. In reality, a portion of it is laminar from the leading edge to 
some downstream position. For this purpose, it was suggested to use 


0.455 A 

(log Re) 2 58 Re 


(10.66a) 


where A has various values depending on the value of Reynolds number at which 
the transition takes place. If the trasition is assumed to take place around a Reynolds 
number of 5x 10 5 , the average skin friction correlation of Schlichling can be written 
as 


C / = 


0.455 


(log Re) 


2.58 


1700 

Re 


(10.66b) 


All that we have presented so far, are valid for a smooth plate. Schlichting used a 
logarithmic expression for turbulent flow over a rough surface and derived 


Cf= 


1.89 


X-2.5 

-1.62 log — 

£ P J 


(10.67) 
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Summary 










Turbulent motion is an irregular motion of fluid particles in a flow field. 
However, for homogeneous and isotropic turbulence, the flow field can 
be described by time-mean motions and fluctuating components. This is 
called Reynolds decomposition of turbulent flow. 

In a three dimensional flow field, the velocity components and the 
pressure can be expressed in terms of the time-averages and the 
corresponding fluctuations. Substitution of these depenent variables in 
the Navier-Stokes equations for incompressible flow and subsequent 
time averaging yield the governing equations for the turbulent flow. The 
mean velocity components of turbulent flow satisfy the same Navier- 
Stokes equations for laminar flow. However, for the turbulent flow, the 
laminar stresses are increased by additional stresses arising out of the 
fluctuating velocity components. These additional stresses are known as 
apparent stresses of turbulent flow or Reynolds stresses. 

In analogy with the laminar shear stresses, the turbulent shear stresses 
can be expressed in terms of mean velocity gradients and a mixing 
coefficient known as eddy viscosity. The eddy viscosity (v,) can be 


expressed as v t =l 2 


du 


dy 


, where / is known as PrandtTs mixing length. 


For a homogeneous and isotropic turbulence, most widely used value of 
mixing length is given by/ =Xy. In this expression ,y is the distance from 
the wall and X is known as von Karman constant (= 0.4). For high 
Reynolds number, fully developed turbulent duct flows, the velocity 
profdeisgivenby 


— =A 1 In T) + Dy 
u t 

where if is the time-mean velocity at any r](=yu T /v) andw T is the friction 

velocity given by V^V/p. The constant.s/1, and/)] are determined from 
experiments whichare2.5 and 5.5, respectively, for smooth pipes. The 
corresponding friction factor (f) is given by 


the expression 


1 

if 


2.0 log 10 (ReV7)-0.8. 


• However, for pipe roughness of varying degree, the following 
recommendation of Colebroolc and White works well 


1 

if 


1.74- 2.0 log 10 



18.7 

ReV7 


where eJR is pipe roughness. 

F C 

• In the range of Re < 10 , the velocity distribution in a smooth pipe is 
given by 
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JL = 8 _ 74 l ^ 


1/7 


• The friction factor in this regime is given by Blasius as 

/= 0.3164 (Rer 0 25 

• The growth of boundary layer for turbulent flow over a flat plate is 


given by 

— = 0.37 (Re v r 1/S 
x 

• The expression for the average skin friction coefficient on the entire 
plate of length/, has been determined as 


Cj = 0.072 (ReJ 


-1/5 


• This result is found to be in good agreement with the experimental 
results in the range of 5 x 10 5 < Re < 10 7 which is given by 


C f = 0.074 (ReJ 


-1/5 


• For turbulent flow over a rough plate, the average skin friction 
coefficient is given by 

/ S—2.5 


C / = 


1.89 +1.62 log — 


p 7 
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Solved Examples 


Example 10.1 Prove that 

1+772 t + T/2 

t-T/2 t-T/2 

where 0 is a constinuous function of E, 

Solution 




t + T/2 


t + T/2 


j 


t-T/2 


t-T/2 


t+T/2 




t-T/2 


dt 
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Changing the order of integration, we can write 

t+T/2 t+T/2 

T 


or 


+1ii t+i // 

J«- I -T 


t-T/2 


t+T/2 


t-T/2 

t+T/2 


t+T/2 

|* 0 d t 

t-T/2 




d<^= d£ 


t-T/2 


t-T/2 


Example 10.2 The well known scientist Theodore von Karman suggested the 
Au! Ay 


mixing length to be l = % 


dru! Ay 2 


. Using this relation drive the velocity profile near the 


wall of a flat-plate boundary layer flow. 
Solution We know 


T, = n t where p, =p/ 2 
dy 


du 


dy 


So, 


T,= pT 


du 


dy 


Substituting von Karman’s suggestion, we get 

__ p x 2 (d« /dy ) 2 (dw / Ay) 1 _ p x 1 {Ad I Ay) 4 


(d 2 w/dy 2 ) 2 


(d ? w/dy 2 ) 2 


An 

4 = 1l J_ 

f d 2 w ^ 

UtJ 

p 2 2 

W) 


Assuming T t = T w and considering u T = ~Jx w l p 

\2 


Au 


" d 2 w^] 

Ul’J 

X 2 

Idv 2 J 


Taking the square root, and applying physical argument that (Au lAy) cannot be 
imaginary, we obtain 


Let 

then 


f Au "I 

2 

U T 

{ d 2 ^ 


- ± -L 


U.vJ 

X 

Uv 2 J 


m = Au/Ay 
AmlAy = ± — m 2 


1 V „ 

-± — y+ C 1 

111 U T 


Integration yields, 
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Using m => °° as v => 0, we get Cj = 0 


Then, 


du 
dy 

Integrating, we obtain, 

At some value 
Invoking this, 


xy 


du n 
since — > 0 
dy 


Thus, 


u = — In y + C 2 
X 

y = y 0 , m = o 

«T 
* 

— = — ln(y-yo) 

«r * 


C 2 = In y 0 


Let us substitute v 0 = ft — order of which is same as viscous sublayer and ft is an 
u T 

arbitrary constant. 

We shall get, thus, 

— = -fln^-lnjS 


or 


— = A 1 In rj + D x 
u T 


This is the universal velocity profile. 

Example 10.3 Using Karman's relation l = X 


du/dv 


d 2 n/dy 2 


, show that the universal 


velocity distribution in a fully developed channel flow (Fig. 10.11) is given by 


^'max u 


In 




where, 2 h is the height of the channel, y is the distance measured from the centre line of 
the channel and X is an empirical constant. The pressure gradient in flow direction is - 
(dp/dx). 

Solution From Reynolds equation, we get 
<?(t,) _ dp 


dy 


or 


r,= 


dx 

dp 2 

y dX; 


y+C, 


At y = 0, T t = 0, that makes Q = 0 


r d]^ 

ydXj 


Thus, 
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and 


T, = 


ydxj 


we get 


y 


h 

From Karman ’s relation, we can write 


= PX~ (du/dv) 
(d 2 u/dy 2 ') 2 


Then 


hpx’ (du/dv) 

To — 


’[d 2 u/dy 2 ^) 


2 = hx 1 (du/dv) 
Uj y(d 2 u/dy 2 f 


Thus, 


d 2 u = ± X_ fh 

dy 2 u T \v 


du 


dy 


Substituting for m = and integrating, 

dv 


-=±2^ Jh>+C 2 


at v = h, m and C 2 = ± 2 — h 

u T 


u T 

X 

U T 


du 

Now, we know that — < 0 for y > 0 and we write 

dy 


du = 


U T f _ 

2 xh J 1 -y[y/h 


dy 


Substituting for 


£,= 1 - 1 j- and integrating, 


at y = 0, u = U m 


d = — [in £ — if;] + C 3 

X 


_ U T 

u = — 


X 


( 


In 


+ Ca 


r = tj + zi_ 

u 3 u max ^ 

X 


Finally, 


^rnax U _ \ 


In 


i-M 
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Example 10.4 During flow over a flat plate the laminar boundary layer undergoes a 
transition to turbulent boundary layer as the flow proceeds in the downstream. It isobserved 
that a parabolic laminar profile is finally changed into a 1 /7th power law velocity profile in 
the turbulent regime. Find out the ratio of turbulent and laminar boundary layers if the 
momentum flux within the boundary layer remains constant. 


Solution Assume width of the boundary layers be a. Then momentum flux is 

x 2 
LL 


where 77 = y/8 

For laminar flow, = 2 n - rr 

U 


H = J u p u a dy = p Ul a 8 J^-^- 


d?7 


i 

A i*m = p^a<5 lam J(4?7 2 -4 tj 3 +t? 4 )dr? 


0 


= put a 8 


lam 


4 3 4 77 

—77 -77 + — 

3 5 


o 2 c 

— yjPtA*> fl ^lam 

For l/7th power law turbulent profile, 

It 


[A 


11 


1/7 


1 

Aurb= pf/ifl<5turbJ(V /7 ) 2 d77 
0 

= puiad^jti 211 dti 


0 


= pUl a8^ 


Comparing the momentum fluxes, 

Slurb - 72 

Slam 105 


^ 9/7 

- 77 

9 


1 q 

= — p Ul a 8 
0 9 P 


turb 


It is to be noted that generally turbulent boundary layer grows faster than the laminar 
boundary layer when a completely turbulent flow is considered from the leading edge. 
However the present result is valid at transition for a constant momentum flow. 


Ex amp l e 10 .5 Air (p = 1.23 kg/m 3 and V = 1.5 x 10 5 m 2 /s) is flowing over a flat 
plate. The free stream speed is 15 m/s. At a distance of 1 m from the leading edge,calculate 
8 andr K , for (a) completely laminar flow, and (b) completely turbulent flow for a l/7th power 
law velocity profile. 
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Solution Applying the results developed in Chapters 9 and 10, we can write for 
parabolic velocity profile (laminar flow) 


8 5.48 du 

- = -S= and 
x ^Rc x dy 


;=0 


Re v = — = 15x1 _ 5 = 1.0 x 10 6 


v 1.5x10 


5.48 

5 = = ■ x 1 m =5.48 mm 


Vl.OxlO 6 


du 

dy 


y =0 


d r?n n 2 l 
5 dv l v 1,7=0 




2 x 1.23 x 1.5 x 10 x 15 
0.00548 


= 0.101 N/m 2 


For turbulent flow, 


or 


0.37 


vl/5 


x (Re A .) 

. 0.370 

o = 


(1.0 x 10 6 ) 1/5 


(from Eq. 10.62) 


x 1 m = 23.34 mm 


5/x = 0.0233 


T w = 0.0225 pUl 


s 1/4 




(fromEq. 10.60) 


T w = 0.0225 x 1.23 x(15) 2 


T w = 0.0225 x 1.23 x(15) 2 


= 0.502 N/m 


A 


1/4 


V X 


\U„x 8 j 

1 


x 


1 


1.0 x 10 0.0233 


1/4 


Turbulent boundary layer has a larger shear stress than the laminar boundary layer. 


Exercises 


10.1 Only write down the option (true/false) or the choice (a, b, c or d) or the 
appropriate conditions. 

(i) For flow through pipes, due to the same pressure gradient, the turbulent 
velocity profile will be more uniform than the laminar velocity profile. 

(True/False) 

(ii) If the mean velocity has a gradient, the turbulence is called isotropic. 

(True/False) 

(iii) + d V = 0 for a turbulent flow signifies 
dx dy 




























The McGraw-Hill Companies 


Turbulent Flow | f~431 | 

(a) conservation bulk momentum transport 

(b) increase in u in positive x direction will be followed by increase in v' 
in negative y direction 

(c) turbulence is anisotropic 

(d) turbulence is isotropic 

(vi) In a turbulent pipe flow the initiation of turbulence is usually observed at a 
Reynolds number (based on pipe diameter) of 
(a) 3.5 x 10 5 (b) 2 x 10 6 

(c) between 2000 and 2700 (d) 5000 

(v) A turbulent boundary is thought to be comprising of laminar sublayer, a 
buffer layer and a turbulent zone. The velocity profile outside the laminar 
sublayer is described by a 

(a) parabolic profile (b) cubic profile 

(c) linear profile (d) logarithmic profile 

(vi) A laminar boundary layer is less likely to separate than a turbulent boundary 

layer. (True/False) 

10.2 Show that, with the help of both the mixing length hypothesis due to Prandtl and 
mixing length law due to Karman (given in Example 10.2), the universal velocity 
profile near the wall in case of a fully developed turbulent flow through a circular 
pipe can be expressed as 

~u = 1 R 
u T k ^R-r J 

where r is the radius of the pipe and K is a constant. 

10.3 Calculate power required to move a flat plate, 8 m long and 3 m wide in water 
(p = 1000 kg/m 3 , p = 1.02 x 10~ 3 kg/ms) at 8 m/s for the following cases: 

(a) the boundary layer is turbulent over the entire surface of the plate 

(b) the transition takes place at Re = 5 x 10 5 . 

A ns. (a) 12.536 x 10 3 W (b) 12.518 x 10 3 W 

10.4 The transition Reynolds number in a pipe flow based on t/ av is approximately 
2300. How does this value can be extrapolated for the flow over a flat-plate if U„ 
in the flat-plate case is analogous to t/ max in the pipe and 8 is analogous to pipe 
radius R7 

Ans. (Re x = 2.116 x 10 5 ) 

10.5 A plate 50 cm long and 2.5 m wide moves in water at a speed of 15 m/s. Estimate 
its drag if the transition takes place at Re = 5 x 10 5 for (a) a smooth wall, and (b) 
a rough wall, e p = 0.1 mm. For water, p = 1000 kg/m 3 and p = 1.02 x 10~ 3 kg/ms. 

Ans. ((a) 411.405 W (b) 806.168 W) 

10.6 In turbulent flat-plate flow, the wall shear stress is given by the formula 
T w = 0.0225 pul 

Two important equations concerning l/7th power law velocity profiles are 

C ft = 2 —— and <5** = —— 8 

* dx 72 

From the above three equations, find the final expression for skin friction 
coefficient (C^). 

Ans. ( C fx = 0.0576 (ReJ^ 175 ) 


v 
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10.7 Water flows at a rate of 0.05 m 3 /s in a 20 cm diameter cast iron pipe 
(s p /D = 0.0007). What is the head (pressure) loss per kilometer of pipe? For 
water, p = 1000 kg/m 3 , v= 1.0 x 10~ 6 m 2 /s. Use Moody’s Chart. 

Ans. (12.2 m) 

10.8 Air (p = 1.2 kg/m 3 and v= 1.5 x 10~ 5 nr/s) flows at a rate of 2.5 m 3 /s in a 30 cm 
x 60 cm steel rectangular duct (e p = 4.6 x 10~ 5 m). What is the pressure drop per 
50 m of the duct? Use Moody’s chart. 

Ans. (217 pa) 

10.9 Water is being transported through a rough pipe line ( u T e p /v= 100), 1 km long 
with maximum velocity of 4 m/s. If the Reynolds number is 1.5 x 10 6 , find out 
the diameter of the pipe and power required to maintain the flow. For water, p = 
1000 kg/m 3 , v— l.Ox 10~ 6 m 2 /s. 

Ans. (£> = 0.454 m,P= 134.113 kW) 

10.10 Modify the friction drag coefficient given by Eq. (10.64) as Cf = 0.074 (Re^) -175 
- Al Re z . Let the flow be laminar up to a distance X ct from the leading edge and 
turbulent for X CI <x <L. Consider the transition to occurs at Re v = 5 x 10 5 . 

Ans. (A = 1700) 

10.11 Air flows over a smooth flat plate at a velocity of 4.4 m/s. The density of air is 
1.029 kg/m 3 and the kinematic viscosity is 1.35 x 10~ 5 m 2 /s. The length of the 
plate is 12 m in the direction of flow. Calculate (a) the boundary layer thickness 
at 16 cm and 12 m respectively, from the leading edge and (b) the drag coefficient 
for the entire plate surface (one side) considering turbulent flow. 

Ans. ((a) 3.5 x 10~ 3 m, and 0.0207 m (b) C f = 3.554 x 10~ 3 ) 

u 

10.12 The velocity distribution for a laminar boundary layer flow is given by — = sin 

u e 


K y 
2 8 


. The velocity at y = k is given by u k . It is assumed that the small 


roughness of height k will not generate eddies to disturb the boundary layer if 
u k k 

- is less than about 5.0. Show that at a distance x from the leading edge, the 

v 

maximum permissible roughness height for the boundary layer to remain 

k A ' xl/4 

undistributed is given by — =- xta 

c (Re) 374 

length of the plate and^4 is a constant. 


— where, Re = , c is the total 







